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Gimballess inertial reference equipment might provide six possible 
advantages over conventional stable member systems, These are savings 
in power, Size, weight, and expense, flexibility of packaging, and ease 
of maintenance. The theoretical problems involved in gimballess equip- 
ment are the question of what inertial instruments measure, and how 
Best to use the information they generate. The engineering problems 
concern the difficulty of torquing gyroscopes with accuracies approaching 
one part in a hundred million, and the design of inertial instruments *o 
withstand the high angular rates they will experience in this application. 
The inertial package to be considered contains three delta-modulated inte- 
grating gyros and three delta-modulated integrating accelerometers, 

Delta modulation is a pulse modulation system which delivers 
Segui pulses, representing increments of the inputeor one of its 
integrals or derivatives, in synchronism with a clock. The analysis 
of delta-modulated loops with either a two- or three-level relay is 
a straightforward matter using z-transforms and modified step-by-step 
transient methods, Application of delta modulation to integrating gyros 
and accelerometers leads to instruments whose output pulses represent 
increments of rotation and velocity. 

Transformation of a vector from one Cartesian coordinate frame 
to another may be represented in three essentially different ways: 
through matrix-tensor operations, four-parameter techniques, and 
vector algebra, Comparison of these to determine which is most 
suitable for machine computation, both from the standpoint of genera*ving 
the transformation parameters and of accomplishing the transformation, 
shows that a matrix transformation, with straightforward generation of 
direction cosines, is best. The most attractive choice of machine for 
the computation is the digital differential analyzer, since the angular and 
linear motion information is in incremental form, High-order integration 
rules with larger angular increments cannot be used to reduce computer 
speed while maintaining accuracy because of an irreducible uncertainty 
in the digital angular information. 
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The transformation computer may be analyzed using frequency domain 
techniques. With these techniques, a signal flow graph representing the 
computer can be developed,which gives great insight into the dominant 
mode of operation of the computer, and which provides an estimate of 
truncation error. Quantization may be treated statistically. 

This investigation has considered one way of instrumenting gimballess 
inertial reference equipment in detail, and has provided a survey of some 
of the possibilities in this field. It has shown the theoretical feasibility 
of such equipment, and has also indicated that the major development 
effort required to permit practical realization of such equipment with 
a specified accuracy must be concerned with inertial instruments designed 
especially for this application. 
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CHAPTER 1 


INTRODUCTION 


1.A Background 

The problem of arriving at a destination has been important to man 
from the beginning. At first, getting there was merely a matter of 
following a line of sight to the destination. Later, landmarks whose 
position with respect to the destination were known were used as 
guideposts in making journeys out of sight of the destination. The stars, 
which can be used out of sight of land, represent an important class of 
these reference points. Celestial bodies can be observed only on clear 
nights, and of coursg, other landmarks, by and large, must be visible 
to be useful. With the invention of radio, and later radar, man-made 
radiation could be used to distinguish landmarks during periods of 
reduced visibility. In modern military and transport applications many 
situations exist in which terrestrial and celestial reference; points 
are either unavailable or unreliable as a basis for guidance. In such 
cases there is a need for navigation systems which do not depend on 
any external information. Inertial navigation systems, which deduce 
the position of a vehicle from on-board measurements, are systems 
of this type. The measured acceleration is integrated twice with respect 
to time to give displacement for an initial position. 

The majority of inertial navigation systems today use two or three 
Ssingle-axis accelerometers mounted on a stable platform to measure 
acceleration. The platform is either kept non-rotating with respect 
to inertial space or precessed so that the input axes of the accelerometers 
are kept coincident with some slowly rotating set of coordinate axes, 


such as local geographic coordinates. The stabilized piatform is 


articulated with respect to its support by a system of gimbals. The 
gyros mounted on the platform detect the angular motion of the platform, 
and the gimbals are torqued in response to gyro signals to drive the 
platform back to its null position. 

Another configuration for inertial navigation systems has the inertial 
instruments fixed to the vehicle and hence is gimballess. Single axis 
accelerometers measure vehicle acceleration and resolve it into vehicle 
coordinates. The guidance equations may be solved in body coordinates, 
in which case gyro information about vehicle rotations is used directly 
in the guidance equations; or they may be solved in an inertial coordinate 
frame, in which case the gyro information is used to generate a coordinate 
transformation which operates on the accelerometer data to give acceleration 
resolved into inertial axes. 

There are six possible advantages available in gimballess guidance 
systems. With the removal of gimbals, weight, spsce, and power 
required are reduced. Further, the not inconsiderable expense of 
gimbals is eliminated. Since the inertial package in a gimballed system 
must assume practically a spherical shape, it is necessary to find a 
spherical space in the vehicle in which to put the guidance system. 

Without gimbals to dictate packaging, much more flexibility in packaging 
is available. Finally,assembly and maintainence of the inertial package 
is easier without the encumbrance of gimbals. 

l.B The Problem of Gimballess Inertial Reference Equipment 

1,.B.1 Theoretical Aspects 

The theoretical problems of gimballess inertial navigation are 
two: the measurement of acceleration in body fixed coordinates and 
the determination of the orientation of the body fixed coordiante frame 
with respect to the inertially fixed coordinate frame from gyro information. 

The three single axis accelerometers measure inertial acceleration 
and resolve it into components along body-fixed axes. Coriolis effects 
are not inherently present in the measured acceleration components; 
they arise only in the intepretation of these components in non-inertial 
reference frames. If such frames are used in the guidance equations, 
the Coriolis terms must be computed separately using angular information 


from the gyros. If the guidance equations are set in inertial coordinates, 
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the three measured components need only be transformed into the inertial 
coordinate system. 

The generation of a vector transformation from gyro information is 
a Straightforward process. The differential equation governing the 
transformation parameters may be derived and instrumented in any 
convenient way. 

1.B.2 Engineering Aspects 

The engineering problems involved in gimballess applications of 
inertial guidance are many, but the two most important involve the 
environment to which the inertial instruments will be subjected and the 
speed of computation required to keep the computed orientation of the 
body -fixed coordinate frame close to its actual orientation. 

While inertial instruments mounted on a stable platform are isolated 
from the angular motions of the vehicle carrying them, they do experience 
all the linear motion and most of the vibrations of the vehicle. Body- 
mounted instruments, in contrast, experience both the angular and linear 
motions, but they may be effectively shock mounted to minimize the 
vibrations present. Although much effort has been expended in developing 
instruments whose performance does not deteriorate under acceleration, 
practically none has been devoted to designs which might experience large 
angular velocities. For this reason no instruments are available at 
present which can be used in a gimballess application. 

In gimballess inertial reference equipment, the gyros are used to 
measure angular orientation of the body-fixed coordinate frame with 
respect to some reference. This may be done either by using a pick-off 
arrangement with free gyros, or by torquing single-degree-of-freedom 
gyros to make them follow body rotations. If free gyros are used, the 
pick:-offs must be accurate to about five thousandths of a degree. The 
required accuracy of torquing for single-degree-of-freedom gyros may 
be shown by the following considerations. If the gyro is to have a drift 
rate of one meru (about 7-10 © radians per second), the uncertainty in 
applied torque must not exceed this ievel. The maximum anguiar velocity 
which a typical vehicle might experience is about one radian per second, 


so that the maximum torque available must be able to precess the gyro 
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at that rate. Comparison of these two numbers shows that for gimballess 
applications, and the above specifications, the torquing must be accurate 
to about seven parts in a hundred million. For higher acceptable drift 
rates or lower maximum velocities, this accuracy requirement is less 
stringent. 

The orientation of the body-fixed coordinate frame must be known 
to a high degree of accuracy (about five thousandths of a degree); thus 
at reasonable angular velocities very high computer Speeds are required 
to keep up with the vehicle rotation. If a single adder is to be used, 
add times on the order of five microseconds per word are required. 
This is near the limit of available performance. 
1.C Scope of the Investigation 

A gimballed inertial reference system delivers acceleration information 
in stabilized inertial coordinates. This investigation will consider an 
inertial reference system with single axis integrating accelerometers 
and single-degree-of-freedom integrating gyros fixed to a vehicle. All 
the instruments will be pulse torqued. The gyro information is to be 
used to generate a coordinate transformation which will operate on 
accelerometer information to give inertially stabilized data. Thus it 
is proposed that the gimballed inertial package be replaced by a body- 
mounted inertial package and a transformation computer. 

This investigation assumes that the guidance equations will be 
solved in inertial coordinates; no examination is made of the possibility 
of doing all or part of the computation in body-fixed coordinates. Further, 
no consideration is taken of the feedback through the guidance computer 
and vehicle autopilot. | 

There are no studies available which compare the attractiveness 
of various coordinate frames available for computation in gimballess 
guidance. For this reason, the present investigation of the performance 
and design of a system which provides inertially stabilized information’ 
from a gimballess inertial package is a useful one. Further, since the 
action of the inertial package and transformation computer is much 
faster than the feedback present, the equipment may profitably be examined 


aS an open loop device. 
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Since the inertial package consists of three delta-modulated singie- 
axis integrating accelerometers and three delta-modulated single-degree- 
of-freedom integrating gyros, the theory of delta modulation and its 
application to inertial instruments wiil be considered in some detail. 
Methods for determining the zero-input behavior of deita-moduiated ioops 
will be reviewed; and the response of such :ioops to impulses and steps 
wiil be found. Criteria for the choice of loop parameters wiil be set 
forth. 

The available methods of transforming a vector from one Cartesian 
coordinate frame to another will be examined. There are essentially 
three ways of viewing the transformation: the matrix of direction cosines, 
the four parameter representations, and the vector representation. The 
properties of each will be considered and the differential equations 
governing the transformation parameters will be derived. Finally, the 
transformations will be compared to determine their suitability for 
machine computation, considering the form of instrument output. 

_ The transformation computer wiil be designed to instrument the 
transformation found # be best for machine computation. The possibility 
of accumulating a number of increments from the gyros and instrumenting 
more complicated algorithms to permit the use of a computer with a 
slower clock rate and reduction of error will be examined. With the 
design of the computer fixed, a detailed analysis will be carried out 
to determine the dominant mode of the computer and estimates of 
truncation and quantization errors. 

Finally a few typical values will be used to illustrate the method 
of determining the basic system parameters and at the same time to 
provide a feei for the order of magnitude of the numbers involved. 


1. D Previous Results 





Deita: modulation was originated as a transmission system (Refs. 
2, 5, 6, 14) with iittle attention being paid to the feedback aspect of 
the transmission loop. Later, with interest developing in contro- 
systems uSing relays, control theory was appiied to deita-modulated 
loops (Refs. 1, 3, 4, 6-13, 15-19) under the names bang-bang, on-off, 


and contactor control systems. The zero-input behavior of these loops 
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received most of the attention, with phase plane methods playing an 
important part (Refs. 8-11). In 1950 a method for determining the exact 
zero input behavior of a loon of any order by means of difference equations 
was presented by Torng and Meserve (Ref. 19) and the z-transform 
equivalent of this method was suggested by Bergen in his discussion of 
their oaper (Ref. 1). Schwartz elaborated Bergen's method and applied 

it to the determination of the loop's response to step inputs (Ref. 15). 
Beyond this, no work has been done previously in determining the response 
of the loop to step inputs. 

The application of delta modulation to inertial instruments first 
took place in 1957, and development work has continued so that these 
instruments are now a very attractive choice for any inertial system 
which requires instrument outputs in digital form (Refs. 4, 7, 8, 17, 18). 

Coordinate transformations have received only cursory examination 
in connection with gimballess inertial reference equipment. Although 
some reports consider the possibility of using transformations other 
than direction cosines, they do not exemine them in any detail. Of 
course, the literature on coordinate transformations iS available in 
almost any text or reference on classical mechanics or mathematical 
physics (Refs. 30-39). 

The analysis of numerical processes has a fairly long history, 
and most techniques are easily applied to digital computers (Refs. 

41, 44, 46, 48, 51-53). However the application of frequency domain 
methods to the analysis of computers is quite rare (Refs. 49-51, 55, 61), 
and thus most analysis to date is quite cumbersome. No attempt has 

yet been made to apply frequency domain analysis to digital computers 

in which the vairable of integration has not been time. 

A number of private companies have worked in the field of 
zimballess guidance, both on their own initiative and under contract. 
Their studies have mostly provided a brief survey of the possibilities 
in this field, but no detailed analyses. Very little of this work is 


available in the open literature (Ref. 66). 


CHAPTER II 


DELTA -MODULATED LOOPS 


2.A Introduction 

Delta modulation is a puise modulation system patented in France 
in 1946 (Ref. 6) and first described in the literature by de Jager (Ref. 5) 
and others (Ref. 14) in 1952 as a system for speech transmission. In 
1957, a similar system was developed giving digital outputs directly 
from instruments and interest in such instruments has remained active. 
(Refs. 4, 7, 8, 12, 17, 18) 
2.B Description of the Loop 

The input is translated into a sequence of pulses, each having the 
same value, but with either polarity being possible. The pulses are 
obtained by means of the negative-feedback loop shown in Fig. 2-1. One 
or the other of the linear plant functions may be unity. When the loop is 
used for transmission, the linear plant in the forward path is usually 
unity, while in instrumentation applications the feedback plant usually 
contains a zero-order hold. In any case, an integration occurs somewhere 
in the linear plant to assure zero average error (Ref. 3, p. 53). The 
relay may have either two levels or three levels. The mathematical 
model for the sampler is an impulse modulator synchronized with a clock. 

Fach impulse appearing at the output represents a unit change in the 
input or one of its integrals or derivatives. This representation of a 
signal in terms of increments leads to the name'dedta modulation". If 
the integration appears in the forward path, each output pulse represents 
a change in the integral of the input. If the integration appears in the 
feedback path, each pulse represents a change in the input itself. Two 


integrations in the feedback path cause each pulse to represent a change 
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in the slope of the input. 
2.C Analysis for Zero Input 


With no input, a loop with a two-level relay will experience a self- 
excited oscillation or limit cycle, and if the zero level is narrow enough, 
a loop with a three-level relay may also have a limit cycle. From physical 
considerations, it is clear that the periods of these limit cycles are 
integral multiples of the sampling.period. Further, if the linear plant 
contains an integration, the number of positive pulses in a cycle will 
equal the number of negative pulses (Ref. 3, p. 53). 

2.C.1 Moding 

Each limit cycle exhibited by delta-modulated loops produces a 
characteristic sequence of output pulses. These sequences are referred 
to as the modes of the system and are very important in the analysis of 
input-output relationships. The modes of a system are designated by 
the output sequence they produce. For example, if the limit cycle ina 
loop with a two-level relay produces a single positive pulse followed 
by a single negative pulse, the mode designation is |-l. Similarly, a 
mode producing two positive pulses followed by two negative pulses is 
designated 2-2. If the loop contains a three-level relay and has a limit 
cycle producing two positive pulses followed by two negative pulses and 
one zero pulse, the mode is termed 2-0-2-l. The zero for the second 
digit indicates that no zero pulse occurs between the positive and negative 
pulses. Note that the mode designations always begin with the positive 
pulses. 

The concept of moding is applicable only to systems with no input, 
since the limit cycle is disturbed by inputs. The description of outputs 
in response to inputs is discussed in Section 2. D. 

2.C.2 Describing Function Analysis 

In his famous paper(Ref. 58) Kochenburger presented a method of 
steady-state sinusoidal dnalysis applicable to non-linear systems in 
which non-linear elements are represented as linear gains. These 
gains are determined as functions of the amplitude and frequency of 


the inputs to the non-linear element. Chow (Ref. 3) and Russell (Ref. 13) 
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first used the describing function to determine limit cycles in delta- 
modulated loops. Their interest was in adjusting the width of the zero 
level in a three-level relay to eliminate limit cycles, but their analysis 
is applicable to any sampled-data loop containing any non-linearity. 

Their method of analysis is the same as the standard describing 
function analysis applied to continuous time systems with the addition 
of a variable delay due to sampling. This variable delay, 6, arises because 
the output pulses can occur only at sampling instants, even though the 
relay had switched earlier. This delay, which ranges from zero to a 
full sampling period, Ts: is easily understood from an examination of 


Fig. 2-2. The marks on the time axis indicate clock pulses. In each 
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Fig. 2-2 Time Delay due to Sampling 


case, the negative output pulse occurs with the second clock pulse shown 
and 6 ranges from almost an entire sampling period in Fig. 2-2(a) to 
practically zero in Fig. 2-2(c). The maximum phase shift due to sampling, 


Qo» is the product of frequency and sampling period, and is given by 


T 
Ss 


yy | 
@_ =OT_ =n aqr- 27a. (2-1) 
Ss 
where T is the period corresponding to the input angular frequency w. 
Standard describing function analysis (Ref. 63, Ch. 10, pp. 559-580) 


is carried out by graphical solution of the equation 


1 + G(jo)N(A, ) = 0 (2-2) 
G(j), the transfer function for the linear part of the system, is plotted 
1 


on a Nyquist diagram with - NW(A, oy’ the negative of the reciprocal of 


the describing function for the non-linear element. Intersection of the 
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two curves indicates a root of Eq. (2-2) and predicts a limit cycle. In 
the analysis ef delta-modulated loops the phase shift due to sampling 


must be added to the plot; this is most conveniently done by adding 9. 


as a phase lag at frequencies on the G(jw) locus for which a limit cycle 
is possible, that is, at submultiples of the sampling frequency. In this 
application, a log polar plot is most convenient for making the Nyquist 
plot. It is termed a ''moding plot". 1 


In Fig. 2-3, a moding plot for G(jw) = is shown. The 


jo(jor +1) 
negative reciprocal describing function for a two-level relay covering 
the whole negative real axis is shown. As can be seen from this plot, 
limit cycles are predicted for angular frequencies of w/2, w./4, and 
w_/6. At w./8 the sampling phase isat most 7/4, which is not sufficient 
to reach the negative real axis. It is clear that limit cycies will not 
exist at lower frequencies; thus the modes predicted are |-l, 2-2, and 
3-3. Fora three-level relay, the describing function reaches a maximum 
when the amplitude of the relay input is V2'times the zero level half-width, 
Aw Thus the reciprocal of the describing function has a minimum, and 
the zero level half-width may be adjusted so that no limit cycle is 
predicted. An example of this calculation is shown in Appendix B. 
2.C.3 Detailed Analysis 

The describing function is only an approximate tool, albeit a powerful 
one. For this reason attempts were made to achieve exact answers. 
The step-by-step transient method was discarded as being too cumbersome. 
Phase-plane methods received more attention, and although limited to 
second order systems became quite elegant (Refs. 9-11). Torng and 
Meserve (Ref. 19) presented a method for determining the exact output 
of the linear plant of a delta-modulated loop at sampling instants using 
discrete orthogonal functions. Bergen (Ref. 1), in his discussion of 
their paper, outlined a similar method using z-transform techniques. 
His method allows the possibility of determining no-transient initial 
conditions for the limit cycle. By manipulation, the state of the linear 
plant at each sampling instant may also be determined. 


Each method assumes a mode at the relay and applies it as an input 
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to the linear plant. The output of the linear plant resulting from this 

input is examined at the sampling instants to see if the signal has the 
proper polarity to sustain the postulated mode. Thus if the 1-1 mode 

be possible, the linear plant output at sampling instants must be a 
repeating sequence of two numbers: the first positive, the second negative. 
Torng and Meserve do this using discrete orthogonal functions in the time 
domain. Bergen operates in the frequency domain with z-transforms, and 
introduces initial conditions in this way. 

The modes which must be investigated are limited, and may be 
determined by a describing function analysis as carried out above (Sec, 
2.C.2). Torng and Meserve list the modes which must be investigated 
for limit cycles with periods of two, four, six, and eight sampling periods 
(Ref. 19, Table II, p. 303). 

For a loop with a three-level relay, choice of the zero-level width 
so that a describing function analysis predicts no limit cycle may be 
made. However, for fast systems i ig > 1), this choice of zero-level 
width may not eliminate low energy limit cycles. To eliminate these, 

a transient analysis is the easiest. The linear plant output is assumed 
to be just at the edge of the zero level of the relay. Thus one feedback 
pulse drives the linear plant back toward null. For no limit cycle to 
occur, the limit of the linear plant motion must not reach the outer 
edge of the zero level, even for very long times. The result of this 
analysis for any linear plant containing an integration is that the zero 
level width should be chosen greater than the product of the sampling 
period and the steady state velocity of the plant under the action of the 
relay alone. 

2.C.4 Compensation 

As Schwartz (Ref. 15) and Sliwkowski (Ref. 16) point out, in limit 
cycling loops inputs may not be sensed until one fuil period of the 
prevailing mode has elapsed. Another view of this problem is that in 
limit cycling, the loop is storing information; the amplitude of the limit 
cycle, indicated by the moding, is a measure of the information storage. 
This information storage should be kept to a minimum if the system 


is to be used to best advantage. In a loop with a three-level relay the 
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width of the zero level is a measure of stored information, and for slow 
systems Cae <1) may become quite large. For this reason it may be 
desirable to compensate the system so that all but the 1-1 mode are 
eliminated, or so that the zero level of a three-level relay is reduced 
to a reasonable size. There are three methods of compensation which 
suggest themselves: de analog, ac analog, and discrete or digital. 

Compensation is best viewed in terms of the moding plot (Fig. 2-3). 
To eliminate the 2-2 mode, more than 90° phase margin must be provided 
for G(jw) at w./4. Similarly to eliminate the 3-3 mode 60° of phase 
margin must be provided at w/6 and so on. With these conditions met, 
the extended G(jw) locus will not intersect the describing function except 
at the 1-1 mode. Fora three-level relay, the zero level width may be 
adjusted so that the reciprocal describing function reaches down to just 
above the 1-1 level. 

De analog compensation methods are covered in detail in most 


elementary automatic. control textbooks (Ref. 57, Vol. I, Ch. 12; Vol II, 





Ch. 4; Ref. 64, Ch. 6; for example). Adjusting the moding plot to achieve 
the objectives listed above is a straightforward task. Ac analog compensation! 
is also treated in many texts (Ref. 63 Ch. 6; Ref. 64 Ch. 6; for example). 
The desirability of this form of compensation stems from the fact that 
the relay is often a phase-sensitive switch operating on information in 
a supressed-carrier wave; thus, if dc compensation were used, the 
Signal would have to be demodulated. However, the signal frequencies 
used in practice are very close to the carrier frequencies, and de- 
modulation is very difficult under these conditions. For this reason, ac 
compensation is quite attractive. It should be noted that the use of low- 
pass to band-pass transformations do not provide sufficiently accurate 
compensation for this application. 

Discrete or digital compensation has received little attention in the 
literature. The advantage of this form of compensation is that there is 
no problem of drift as compared to analog compensation, and non-linear 
or logical compensation is possible. Most books on sampled-data systems 
include sections on the techniques of linear discrete compensation and 


Tou (Ref. 62, pp. 444 ff) lists three methods for realizing linear discrete 
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transfer functions: computers, delay lines, and pulsed RC circuits. 

The point in the delta-modulated loop at which the signal is discrete 
is after the sampler; it is here that the discrete compensator must be 
inserted. Ifa pulsed RC circuit is used, a second sampler must follow 
NG 

Linear discrete compensation is best viewed in the G(z) plane. The 
locus of the linear plant has extensions similar to those in the moding 
plot, because of sampling delays. The compensation must be chosen so 
that the extended G"(s) locus does not intersect the negative reciprocal 
describing function locus. 

Logical discrete compensation, using branching operations which 
influence the relay on the basis of past output pulses, has been tried only 
once (Ref. 20) with little success. It would seem that there are great 
possibilities for this type of compensation. 

2.D_ Input-Output Relationships 

As explained above, each output pulse represents a unit change in 
the input or one of its integrals or derivatives. Thus with a finite input, 
the algebraic sum of the output pulses should average to the input or one 
of its integrals or derivatives or, in the case of zero input, to zero. 

The system overloads when the rate of change of the input (or one of 

its integrals or derivatives) exceeds that value which requires all the 
output pulses to be of one sign (Ref. 2, p. 145). Therefore a convenient’ 
way of expressing the input is as a fraction of the maximum permissible 
signal. 

Since the position and number of integrators in the loop change input- 
output relationships, a general analysis is very cumbersome. For this 
reason, the analysis discussed here considers a loop with one integrator 
init. The linear plant in the feedback path is moved to the forward path, 
with the exception of zero-order hold, if one is present, which is left 
in the feedback path. This requires the reciprocal of the moved plant 
to be placed in front of the loop. Fig. 2-4 shows the rearranged loop. 
The dotted lines around the hold indicate that it may or may not be 
present. Gi(s) is the part of the linear plant moved out of the feedback 


path. The analysis which follows ignores the 1/G5 (s) in front of the 
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loop since it tends to confuse the issue. 
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Fig. 2-4 Delta-Modulated Loop, Rearranged 
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2.D.1 Transient Response 

The following analyses are exact determinations of the inpulse and 
step responses of delta modulated loops. With the exception of Schwartz's 
work (Ref. 15) on step inputs for loops with two-level relays, they have 
not been previously described. 
2.D.1l.a Impulse Response 

The height of an impulse (infinity) exceeds the capacity of the system; 
the best output makes up the integral of the impulse a step at a time to 
within the quantization error of the loop. Clearly if the value of the 
impulse is too low, no output will appear, except for the limit cycle. 

The analysis of the loop's impulse response is best done by a step- 
by-step transient method. The applied impulse drives the linear plant 
Some distance from null; calculation of this excursion in the presence 
of the feedback signal is a simple task. The linear plant output at 
sampling instants is examined for the time at which the feedback signal 
changes sign. If the loop has a two-level relay, the state of the linear 
plant is compared with the set of initial conditions corresponding to 
each mode to determine into which mode the system will settle. \ Usually 
this will be the highest mode. The output of the linear plant plant ina 
loop containing a three-level relay is examined to See if it will cause 


the relay to have an output of opposite sign from the input. If not, the 
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computation is complete; if so, it must be carriedfurther. 
2.D.1.b Step Response 

A step input of magnitude D, where D is the drive level of the relay, 
will saturate the system. For this reason the magnitudes of step inputs 
are given as a fraction, r, of D. 

For a periodic sequence to appear at the output, the input must be 
expressed as a rational fraction of the saturating input. This is nota 
serious limitation, since any number less than one may be well approximated 
by a rational number with a small denominator 

Suppose the input is a rational fraction, say U/U, of the saturating 
signal. Then the output will be a sequence of NW pulses. For a loop 
with a two-level relay, there will be Nt = (W+U)/2 positive pulses and 
N = N(W-U)/2 negative pulses. This may be shown by a short computation. 
The relationship between U, W, N’, and N is 

N’-N7 _ UL 


N+N7 OW pe2) 
Rearrangement of this ratio yields 

NO. Wty 

N Wee) (2-4) 


A loop with a three-level relay may have N° zero pulses. In this case, 


+ - 
the required relation between U, W, N , N , and N° is 








NENG. OU 
N'+N-+N° W (2-5) 
If U is positive, and there are no negative pulses, then 
~- 
Ne ey 
N° w-U (2-6) 


With ratios, of positive, negative, and zero pulses determined, it only 
remains to determine the organization of these pulses within the sequence. 
This may be done by determining bounds on the number of consecutive 
pulses of one sign under various conditions. Different bounds are of 
interest depending on whether the loop has a two-level or a three-level 


relay. 
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For a loop with a two-level relay, the bounds that are important are 
the maximum number of pulses of one sign that are possible and the 
minimum number of pulses of one sign that are possible following a 
given number of pulses of the opposite sign. 

The upper bound is found by observing the greatest possible excursion 
from null and counting the number of clock periods required to drive the 
linear plant back across null. The maximum excursion from null occurs 
when the linear plant output travels away from null at maximum rate for 
a full sampling period. 

The lower bound is found by assuming minimum initial conditions. 
The minimum displacement is clearly zero plus, and the minimum 
derivatives result from maximum velocity at the previous switching 
time a given number of sampling periods earlier. Under these initial 
conditions, the number of sampling periods required for the system to 
cross null is found. 

The bounds to be determined in the case of a loop with a three-level 
relay are the maximum number of pulses of the same sign as the input 
and the maximum number of opposite sign. 

Determination of the maximum number of pulses of the same sign 
as the input proceeds from the assumption that the linear plant has 
travelled at maximum rate under the influence of the input alone for a 
full sampling period. At this point, positive pulses begin to drive the 
plant back to null. The pulses required to bring the linear plant output 
within the zero level are counted. 

The possibility of pulses of a polarity opposite to the input may be 
examined by a similar analysis. It is assumed that the linear plant is 
driven across the zero level under the influence of the feedback signal 
and input for a full sampling period starting from just outside the zero 
level. The feedback signal switches off, and the system output is 
examined at successive sampling times to see whether it gets outside 
the zero level on the other side. 

Examples of the above computations are shown in Appendices A and B. 


With these bounds set, the sequences having the required average 
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may be written down. These sequences are designated in a way similar 
to that used for modes. For example, a sequence from a loop with a 
two-level relay consisting of three positive pulses followed by a negative 
pulse is designated 3, 1; one having four positive pulses, one negative 
pulse, two positive pulses, and another negative pulse is designated 

4,1, 2, 1. A loop with a three-level relay having an output sequence of 

a positive pulse followed by a zero pulse would have that sequence de- 
signated 1, 1, 0, 0. Nothe that each of these sequences has an average 

of one-half. The permissible sequences, which are few because of the 
elimination process just described, may be tested by Bergen's method 

to see if they can actually occur. The sequence, with the input, is assumed 
to drive the linear plant, and the signs of the linear plant output are 
examined at sampling times to see if the postulated sequence is sustained. 

With inputs smaller than about one-tenth of the saturating signal, 
the period of the output sequence becomes rather too long for the analysis 
above to handle conveniently, and a modified step-by-step transient 
procedure is preferable. 

In a loop with a two-level relay and a small input, the prevailing mode 
will continue for a number of cycles. The feedback signal has zero average 
over a full period of the mode, and its effect may be ignored. This is 
perhaps best visualized in phase space with the mode trajectory moving 
slowly through the space under the influence of the input. When the 
position coordinate of one of the sampling points of the trajectory moves 
to a new region of phase space (from positive to negative, for example), 
the prevailing mode will change. The mode changes either by dropping 
a pulse or by adding a pulse. In any case, the total state of the linear 
plant after the mode changes may be found and compared with the initial 
states of the other modes. Thus the new mode of the system is found. 
The trajectory of the new mode is assumed to move under the influence of 
the input, with a new initial position and an initial velocity which is the 
velocity of the plant relative to the new prevailing mode. ‘This process 
is continued until the entire sequence is determined. 

For the loop with a three-level relay, the appropriate calculation 


is a step-by-step transient response. Assume an initial condition, and 
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follow the plant through until the sequence is determined. In almost 
all cases only one positive output pulse will occur between zero pulses. 
2.D.1.c Effects of Imperfect Integration 

Thus far, it has been assumed that the integrator has been perfect. 

In some cases, the integrator is in reality a first order lag with a long 
time constant. This may be viewed as a small elastic restraint. With 
an elastic restraint present, there is a limit below which inputs cannot 
be detected. This is also best seen in terms of the mode trajectory in 
phase space. If the steady-state response of the linear plant to a small 
input is insufficient to move the trajectory far enough to cause the mode 
to change, no input will be detected. Of course, with a three-level relay, 
this critical value of linear plant output is the edge of the zero level. 
2.D.1.d Ramps and Higher Order Inputs 

Ramps and higher order inputs (e.g. parabolic) cannot be followed 
by the loop over infinite time since the amplitude of such a signal increases 
without limit. A sawtooth wave may be analyzed in terms of the frequency 
response described below. For ramps which rise to a certain value and 
then remain constant, a step-by-step transient computation seems to be 
the most convenient. The bounds on number of pulses are useful in this 
computation. 

2.D.2 Frequency Response 

since the frequency response method of servomechanism design is 
familiar to most engineers, the following modification of the describing 
function is presented. 

A sinusoidal signal is assumed at the relay. The standard describing 
function is used to represent the relay, assuming enough low pass filtering 
in the loop to make the describing function valid. The fundamental of the 
relay is passed through the sampler and hold, giving the feedback signal. 
With the input to the relay, which is the output of the linear plant, known, 
the error signal can be determined. The input to the loop is the sum of 
the error and feedback signals. With the relay output taken as the loop 
output, the describing function for the entire loop is determined as the 
ratio of reiay output to loop input. 


Three comments must be made about the above analysis. First, 
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the output of the loop has been taken as the relay output, with no limit 
cycle present. Thus, this analysis, for frequencies approaching the 
limit cycle frequency, is open to question. Second, there is a minimun 
amplitude of signal which can be detected,corresponding either to the 
limit cycle amplitude or the zero-level width of the relay. Third, the 
bandwidth of the loop is limited by the sampling frequency, and signals 
above this frequency will not be passed; instead they appear as steps. 
(Ref. 19). 
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CHAPTER III 


DELTA-MODULATED INSTRUMENTS 


3.A General Comments 

The delta-modulated loops discussed in Chapter II are to be applied 
to two inertial instruments, the gyroscope and the accelerometer. 
The gyro to be used is the single-degree-of-freedom floated integrating 
gyro described by Draper, Wrigley and Grohe (Ref. 21). The 
accelerometer is a single-axis floated integrating pendulum type 
described by Hutchings and others (Refs. 4, 7, 8, 12, 17). 
3.A.l1 Brief Description of the Instruments 

The integrating gyro, when behaving in an ideal manner,is 


described by the following differential equation: 


JAoa a Oye 


OA Ae Vta id Mie (3-la) 


where J is the moment of inertia of the gyro float about its output 
axis; 
C is the damping coefficient of the fluid damper; 


Hep is the spin angular momentum of the gyro; 
Wia is the angular velocity of the gyro case about its input axis, 


with respect to inertial space; 


Mis is the torque applied to the instrument float by the torque 


generator; 


and Aon is the angle through which the float has turned with 


respect to the case avout the instrument output axis. 


A dot over a quantity (e.g. Aga) denotes time differentiation. 
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Fig. 3-1 Delta-Modulated Instruments 
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Eq. (3-la) may be rewritten as 


° : H M 
TAga + 40q = Wy t (3-1b) 
where 7 = J/C is the characteristic time or time constant of the gyro. 


The differential equation describing the accelerometer is 


enc 7 Ceo yan ae (3-2a) 
where the symbols above are similarly defined; where P is the 
pendulosity of the accelerometer; and ara is the non-gravitational 
acceleration of the instrument with respect to inertial space along 
its input axis. 


Eq. (3-2a) may be rewritten as 


on + Aga 7 1a +e (3-2b) 
where again 7 = J/C. 

The torque generators are fed from a constant current source, with 
the current switched between two opposing windings to produce torque 
of different signs, or switched to a dummy load or off to produce zero 
torque. The required sign of the torque is determined at sampling 
times and is either opposite to the sign of the float angle or zero, 
depending on the logic of the loop. Therefore, the torque generator 
acts as a zero-order hold in the feedback path. 

3.A.2 Form of the Loops 





The delta-modulated instruments may be represented as shown 
in Fig. 3-l. These loops are identical in form to each other and to 
the loops analyzed in Appendices A and B. The output pulses represent 
changes in the integral of the input, so that the gyro's output pulses 
indicate changes in the angular orientation of the gyro case about its 
input axis, and the accelerometer's output pulses indicated changes 
in velocity along the accelerometer's input axis. 
3.B Delta-Modulated Gyros 
3.B.1 Detailed Description 


A detailed physical description of the single-degree-of-freedom 
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floated integrating gyro is presented in Reference 21. Current 
practice is to center the gyro float in the case with a magnetic 
suspension system similar to that described in Reference 26, as 
well as to float it. Nominal values of the characteristics of some 
gyros are given in Reference 25. Fig. 3-2 shows a schematic 
drawing of the instrument. 

The differential equation governing the gyro may be derived 
from Euler's equations of motionfor arigid body about principal 
axes. It is assumed that the bearings of the gyro are infinitely stiff, 
so that only motions about the output axis need be considered. The 
bearing stiffness assumption is quite well satisfied for gyros even 
under high angular velocities. The set (IA, OA, SRA) forms a right- 


handed triad. The equation is 


Map = ~Hep!Wq 005 Aga - Wop, Sin Aggy) + HAoat W..) 
+(I, - Te ee - Wages) sin Aj, COS Aon 
+(la> Igy eos“Agy - sin” Agy) Wy, Wop (3-3a) 
Mi op = Mi ~ en (3-3b) 


where M. is the torque applied to the instrument float about 


Ppp 


its output axis; 


WwW is the angular velocity of the gyro float about its spin 


SRA 
reference axis, with respect to inertial space; 


vo is the angular velocity of the case about its output axis, with 


respect to inertial space; 


La is the moment of inertia of the gyro float about its input axis; 
Ian is the moment of inertia of the gyro float about its spin axis. 


These equations may be rewritten as 
» ep ae + 
TAog * Aga - > Wgpa SM Aga 7-H Wry 608 Aggy 
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Fig. 3-2 Pictorial Diagram of Floated Single-Degree-of-Freedom 


Integrating Gyro 
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The last two terms on the right side of Eq. (3-3c) are anisoinertia 
effects and are negligible in practical gyros. They will not be 
considered further here. 
Since Aoa stays well below a degree in practice, the standard 


small angle approximations may be made in Eq. (3-3c): 


¢o9¢g eo H H e M 
_ _Sp - Sp P tg 


We Oo Be (3-3d) is identical with Eq. (3-lb). 


(3-3d) 


or Weopa = 


3.B.2 Rotational Coupling 

Eq. (3-3d) shows that rotations of the gyros about all axes affect 
its motion. The coupling effects of these motions will be considered 
in detail. 
3.B.2.a Non-Reversing Rotations 

Examination of Eq. (3-3d) shows that an angular velocity about 
the spin reference axis appears as an elastic restraint term in the 
gyro equation. This elastic restraint causes imperfect integration in 
the loop. As noted in Chapter II, Section 2.D.1.c, imperfect integration 
in the loop produces a dead zone, and inputs within this dead zone 
cannot be detected. This is a significant phenomenon in gimballess 
inertial reference equipment, since large angular velocities may 
occur about any axis. 

Angular acceleration of the case about its output axis produces 
a torque on the instrument float of the inertia reaction category. As 
can be seen from Eq. (3-3d), this torque is proportional to the gyro 
time constant and acts as a spurious input to the gyro. 

These torques may be compensated out either directly in the 
gyro or in the computer, if it is found that they seriously affect the 
performance of the gyro, in a manner similar to that used for 


compensation of compliance effects in present gyros. 
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So SA ®) Vibratory Rotations 


Vibratory rotations of the gyro about one axis at a time do not 
affect the performance of the instrument. When this type of motion 
occurs about two axes at a time, however, it can produce a drift. 

Consider sinusoidal motion of the case about the spin reference 
axis and the input axis. With sinusoidal motion, the gyro equation of 


motion, ignoring torque generator torque, is 


_ sp _ sp 
TAon + Aon == Wopnan sin (27 fopatAcoa = Wia sin (27 frat + +) 


(3-4) 
A solution for this equation using series techniques, which coverges 


for Wopan H../C <1, has been given by Haff and Meltzer (Ref. 27, Ch. 5). 


A further discussion of the series solution and its convergence is given 
in Reference 29. 


Haff and Meltzer show that if nfopan = mf a: a constant drift is 


predicted by the series solution which they propose. A similar effect 
occurs if a vibration about the output axis occurs with the vibration about 


the spin reference axis. Further, if foRA is harmonically related to 


the clock frequency, drift may also be present, as can be seen by writing 
the Fourier series expansion for the torque generator torque. 

Angular vibrations of the gyro about any two axes produce a conical 
motion. Computation of the orientation of the gyro may be made in such 
a way as to introduce a constant angular velocity about the cone axis in 
the computed orientation, when none actually exists. Haff and Meltzer 
present a general treatment of this effect for a stabilized platform 
(Ref. 27, Ch. 3, Sec. 3.4). It is clear from their discussion how the 
computation implied by a stable platform leads to this error. This 
coning error does not exist in gimballess systems. 

3.B.3 Linear Motions 

Many writers (Ref. 22, for example) have treated the effects of. 
linear motions on gyros and they will not be considered here. 
3.B.4 Choice of Loop Parameters 


While the parameters of the loop must be chosen so that production 
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gyros and available electronic components may be used, a certain 
amount of latitude exists in choosing some parameters. 

The spin angular momentum of the gyro is chosen on the basis 
of the allowable drift rate. The higher the drift acceptable, the 
lower the angular momentum may be. 

The torque level, M, required in the torquer depends on the 


maximum angular velocity that the system must follow, Wis eess The 


product of angular momentum and W gives the lower bound on M. 


max 
For gyros with drift characteristics acceptable in inertial guidance 
applications, this bound may be quite high, and torquers other than 
microsyns may be required. 

The clock frequency may be determined from the maximum 
angular velocity expected and the increment size in the loop, AQ@,. 


The quotient W_., /A4@ gives the lower bound on clock frequency. 


If an ac signal generator is used, the carrier must be a multiple of 

the clock frequency. This requirement follows from the desirability of 
peak detection in the switch. If the carrier frequency is not a-multiple 
of the clock frequency, the phasing of the clock and the carrier will 
not remain constant, and the clock pulses cannot be adjusted to appear 
at the peaks of the carrier sinusoid. 

Since it is desirable to have the uncompensated mode of the loop 
as low as possible, or alternatively to have the relay zero-level width 
as narrow as possible, the gyro time constant, J/C, should be as 
low as possible. Since the polar moment of inertia of the float is 
usually fixed by the choice of angular momentum, the time constant 
depends only on the damping coefficient, C, which can be varied by 
varying the operating temperature of the instrument and the clearances 
between float and case. The higher the damping, the lower the time 
constant, but the smaller the angle through which the float travels. 
There is a certain value, Pak, below which the polarity of the signal 
cannot be determined, so that an indefinite increase in damping is not 


practicable. A short computation, suggested by J. E. Miller of the 
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M.I. T. Instrumentation Laboratory, gives an expression for the 
optimum damping in a loop with a two-level relay, considering only 
the dynamics of the gyro, which usually dominate. 


The differential equation governing the gyro is 


> png * 


er OL ap 1A te (esa) 


Integration of this expression over an integral number of sampling 


+ CA 


periodsz and assumption of zero initial conditions for Ao A gives 


We CA =n (6s. -@ 


OA ony Hgp?14 ° Find) (3-5b) 


where aN is the angle through which the gyro case has turned about 


its input axis; and 


0. q is the indicated value of OTA? obtained by counting output 


pulses. 


The quantity (954 - 9...) may be recognized as information stored 


ind 
in the float (see Ch. II, Sec. 2.C.4). When the torque switches sign, 
: ' : n 
the value of Ao a is Ea and the value of Ao A is M(1-a' lta re. 
T 
—s 
Here @a=e 7 andthe loop is assumed to be moding n-n. Thus at 


switching, 


_ JM(-2') ,C_ 
n A. 
Hype ate Hsp a (3-5c) 
By differentiating Eq. (3-5c) with respect to C, and setting the 


Ss 


derivative equal to zero, an extreme value of I. may be found, which 


can be shown to be a minimum. The dependence of a and non C may 
be neglected. 
min * 9 


Sp 


dl. . JM(l-e"), A 
S = weer 
qc H Cite") F 
sp 





or 
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C = l(1-a") JM 
n 
To show that this value of C indeed provides a minimum in I. 


d°I 


S is computed: 





dC 


n 
dT, 259 JIM(l-a), 4 


dc H C2(1+e") 
sp 


which indicates a minimum. 
This optimum value of C has been chosen to minimize stored 


information. The corresponding time constant is 


n 
J (lta) J a 


C (l-a )M (3-7) 
For loops with a three-level relay, the damping should be chosen 
so that 
a Amin 8) 


The expressions for A. are derived in Appendix B, Sec. B.1. 


It is possible to show a relationship between the elastic restraint 
present in a gyro and the threshold of the instrument. The differential 
equation of the gyro is 

JA aa OF = + KA 


=H W,, +M 
sp t 


OA OA OA TA g (3-9) 
where K is the elastic restraint acting on the float from all causes, 
including electromagnetic reaction torque and the effect of angular 
velocity about the spin reference axis. If the inertia-reaction torque 
is neglected (it is usually much smaller than the damping torque) and 
the average angle of the float, ae is considered, limit cycles may be 
ignored (See Ch. II, Sec. 2.D.1.b) Under these assumptions, the solution 
of Eq. (3-9) is 

_K 


H_ W 
AA, = sp _ih (cule) 


OA (3-10) 
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In time T, the average angle of the float reaches a value A Ag? which 


represents one increment in the loop's output. This value is given by 


_M 
Ang = cls 


Writing Eq. (3-10) for t = T and solving for T yields 


Cc AQ 
T=-—~0n (1-7-—y— K) 
- SD IA (3-11) 


T becomes very large as (A 9 / Bs Wha)K approaches unity. For this 


term greater than or equal to one, Tis not defined. Physically, this 
corresponds to the case where elastic restraint completely masks 
the input (See Ch. II, Sec. 2.D.1.c). Using only the first and second 
order terms in the series expansion of the logarithm, for small 


K, gives 


C9 aes 
eS K) 
H = W 7 ad's Teas), 
7 as: Ques Gay 
Since MT, = Hp O? this may be written 
p = Ag (145 are +S) 
IA IA (3-13) 
Thus the percentage error in time of occurrence of an output pulse 
is 
E€ = é (aul a nS 
a (3-14) 


By manipulating Eq. (3-14), many interesting relations between the 
time error and the loop parameters can be seen. The error in time 
decreases with elastic restraint and maximum input angular velocity. 
It is inversely proportiqnal to the actual input angular velocity and 
the square of the damping coefficient. Since the optimum damping 


coefficient depends inversely on Ain? decreasing AL * decreases 


i 


er 


ol 


3.C Delta-Modulated Accelerometers 
3.C.1 Detailed Description 

A detailed physical description of a pulse-restrained accelerometer 
is given in Reference 8. Present practice is to center the float in 
the case by a magnetic suspension system described in Reference 25. 
Reference 4 gives some parameter values for a typical accelerometer. 
Fig. 3-3 shows a schematic drawing of the instrument. 


The accelerometer's differential equation is 


JAoat CA + Pal = Pa 


cos A - JW 


OA PRA 82 Aga 1A OA i 


Joule + M, 


1a ~ tpa)WyaW pa g (3-15a) 
where PA and PRA refer to the pendulum axis and pendulum reference 
axis respectively. (IA, PRA, OA) is a right-handed triad. 

The last term on the right hand side of Eq. (3-15a) is an anisoinertia 
term which is negligible in practical accelerometers and it will not 
be considered further here. Making the small angle approximation 


for Ao jn and dividing by C gives 


eae : - iP : 
Ton *Aoa-G @pratoa =e *1a ~ oat Ce 
(3-15b) 
When 2.55, = W,, = 0, Eq. (3-15b) becomes Eq. (3-2b). As in the 
case of the gyro, infinite bearing stiffness has been assumed, and 
torques considered only about the output axis. 
3.C.2 Rotational Coupling Effects 
As can be seen from Eq. (3-15b), only rotations of the instrument 
about its output axis affect it, because of the assumed stiffness of 
the bearings. As in the case of the gyro, angular accelerations of 
the case about the output axis look like inputs. 
If the case motion about its output axis is vibratory and a vibratory 
linear motion along the pendulum reference axis occurs, an equation 
of the form of Eq. (3-4) governs the accelerometer motion. As noted 
above, if the frequencies of the two motions are harmonically related, 


a drift in the loop's output may appear. 
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3.C.3 Linear Coupling Effects 

Eq. (3-15b) shows that an acceleration along the pendulum reference 
axis appears as an elastic restraint. Imperfect integration results 
and a dead zone appears in the instrument. This is important in any 
application of this accelerometer since it may experience acceleration 
in any direction. 

If a vibration occurs along the pendulum reference axis in 
conjunction with a vibratory input, a drift may result as described 
above. A similar situation may arise if vibrations along the pendulum 
reference axis are related to the clock frequency. 

3.C.4 Choice cf Loop Parameters 

Parameters for delta-modulated accelerometers may be chosen 
in a way Similar to that used when choosing gyro parameters. Those 
which may be chosen are pendulosity, torque level, clock frequency, 
and damping. 

The pendulosity of the accelerometer is chosen on the basis of 
the level of torque uncertainties expected. The lower the torque 
uncertainty, the lower the pendulosity may be. This choice involves 
the same sort of compromise as the choice of angular momentum for 
a gyro. 

The required torque level is the product of the pendulosity and 
the highest level of input acceleration expected. This value is practically 
always lower than the corresponding value for a gyro. 

The clock rate is determined from the maximum acceleration 
expected and the velocity increments used. Their quotient gives the 
lower bound on clock frequency. 

The choice of damping coefficient for the accelerometer is made 
just as in the case of a gyro. Fora limit cycling accelerometer, the 


optimum damping is 


SS 
Nn 
., /{i-@ JJM 


C 0 
\ GQt+a)A.. 
min 


(3-16) 


which gives a time constant 


o4 









(Ite)JA_. 
min 


(1-o")M (3-17) 
For an accelerometer with a three-level relay, the damping coefficient 
should be chosen so that 
Ay = Amin (3-18) 
The remarks about elastic restraint made with respect to the 
gyro apply to the accelerometer as well. The equation corresponding 
to Eq. (3-14) is 


cn ob AY 


T a 


QA 


(3-19) 
3.D Non-Ideal Components 

Up to this point, very little has been said about the non-ideal 
behavior of any components. The following sections consider some of 
the more serious hardware problems. 
3.D.1 Gyros 

In most applications, floated, single-degree-of-freedom integrating 
gyros do not experience significant angular velocities about any axis. 
In any environment in which they do, their performance will deteriorate, 
simply because they were not designed to take such punishment. A high 
angular velocity environment affects present-day gyros in two ways. 
First, they are not structurally strong enough to withstand the bending 
and torsion they experience. When torque is applied by the torque 
gererator, the float experiences a torsional and bending stress because 
of the high spin angular momentum it has. In gimballess applications, 
the high torques involved may literally tear the instrument apart. 
A second effect of high angular velocities is high bearing loading. The 
bearings must be able to supply reaction forces strong enough to counter 
a torque equal to the torque generator torque. This cannot be done for 
long periods of time, but in a floated instrument the floatation fluid 
provides temporary additional’ bearing support, so that for short times 
the float will stay réasonably well centered. While at the present time 


floatation provides necessary additional support for the instrument, 
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magnetic supports are becoming capable of handling the entire load 
alone. 
3.D.2 Torque Generator and Current Source 

In the analysis above, it has been assumed that the time constant 
of the torque generator was negligibly small. This assumption is 
violated as the torque from a given instrument is increased. Torque 
rises as the square of the input current. To increase this current, 
either the voltage must be increased or the impedance of the torquer 
decreased. Increasing the voltage is not an attractive solution, since, 
in general transistors drive the torquer, and their voltage limit is 
quite low. Lowering the impedance requires lowering the resistance, 
which in turn raises the time constant of the torquer, L/R. As the 
time constant of the torque motor increases, the effect of the lag 
makes itself felt in the moding plot. The value of the torquer time 
constant above which the lag must be considered is about The -$T.. 
More torque can, of course, be obtained by increasing the size of the 


torquer, but it is desirable to have small instruments. 


Microsyn torquers (Ref. 28) have been successfully used in pulse 
torquing applications. However, they are not high torque instruments, 
and their dependence on current squared makes them non-linear when 
used with a three-level relay. 

Flapper-type torquers (Ref. 23) provide about ten times as much 
torque aS a corresponding size microsyn. Their main disadvantage 
is the presence of a fairly large negative elastic restraint, which 
makes the torquer positionally unstable. This might be eliminated by 
providing a positive mechanical or electrical spring, or using ac 
voltage for torquing and tuning out the restraint term. The ac torquer 
provides a double frequency term as well as a constant average, so 
that the frequency used must be well above the clock frequency, and 
above the mechanical break point of the system. 

Moving-coil, or d'Arsonval, instruments also provide large torques 
for small size. This type of torquer usually uses a permanent magnet 
field. The difficulty with permanent magnets is that measurement 
techniques are not yet accurate enough to determine the stability of the 
field. 


A problem that all three kinds of torquers share is the stability 
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of their sensitivities. Good operation of the instruments requires that 
the same angular impulse be applied to the float for each electrical 
impulse applied to the torquer. If this is not the case, the instrument 
scale factor will vary with time, and errors will appear from this 
source. 

A related problem is the variation in torquer time constants with 
direction of torque. If the time constants are different, more impulse 
will be applied in the direction of the smaller time constant. Eventually 
a change in the output will be indicated although the input has been zero. 

To achieve the required accuracy in torquing .‘a: VELP poet oan 
stable current source must be used. If a three-level relay is used, 
either the current must be turned off, or it must be fed into a dummy 
load. Turning off the current, while achieving great.savings in power 
causes Stability problems in the current source above those already 
present. Ac current sources can be made more stable, but ac torquers 
are less sensitive than de torquers. 

3.D.3 Signal Generator 

With ac signal generators, the output is a supressed-carrier wave. 
Thus, careful attention must be paid to the effects of circuitry which 
operates on the signal generator signal so that the gain for each 
sideband frequency is the same, and the phase shifts at each sideband 
are the negative of each other. If this is not the case, a quadrature 
component,which will cause a larger uncertainty in the position of the 
float with respect to the case, will be introduced into the signal. The 
end result will be an increase in the minimum angle of the float which 
can be detected. For this reason, the low-pass to band-pass trans- 
formations commonly used in ac compensation are not sufficiently 
accurate to be used in this application. 

As noted above (Sec. 3.B.4) when ac signal generators are used, 
the carrier and clock pulse phases are adjusted so that the clock pulses 
occur at the peaks of the carrier sinusoid. This adjustment is usually 
made by shifting the ac signal at either the primary or the secondary 


winding of the signal generator. For the reasons noted above this 


of 


phase shifting should be done on the primary to avoid adding a 
quadrature component to the signal generator and increasing the 
minimum discernable angle. | 
Flapper-type (Ref. 24) and d'Arsonval signal generators provide 
attractive alternatives to microsyn signal generators, because of 


their small size. 
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CHAPTER IV 


KINEMATICS OF ROTATING COORDINATE FRAMES 


4,.A Introduction 

In a conventional inertial guidance system, the inertial space reference’ 
is preserved physically by the gimballed stable platform. The gimbals 
isolate the platform from the motion of the vehicle through servo drives, 
operating to null the output of gyros mounted on the stable platform. In 
a gimballess system, the inertial space reference may be preserved 
analytically in a transformation computer. This computer operates on 
a vector measured in body-fixed coordinates (the AV vector) to give the 
components of the vector in an inertially fixed frame of reference. The 
computer uses angular information obtained from gyros mounted on the 
vehicle to generate the rotational transformation. 

In this chapter, the properties of the transformation and the relations 
between the angular velocity of the body-fixed coordinate system with 
respect to inertial space and the rates of change of the transformation 
quantities will be adduced. 

There are essentially three methods of representing rotational trans- 
formations. Probably the most familiar is the matrix of direction cosines. 
This matrix, or alternatively, this second-order tensor, may have its 
elements written in terms of the nine direction cosines explicitly, in 
terms of a set of Euler angles, or in terms of the axis and magnitude of 
the rotation. Quaternions, invented by Sir William Hamilton in 1853 (Ref. 31), 
use four parameters to describe a rotation. The vector to be transformed 
is cast in quaternion form, pre-multiplied by a quaternion representing 


a rotation, and post-multiplied by the conjugate quaternion. The result 


is a quaternion representing the rotated vector. A variation on this, the 
Cayley-Klein parameters, wes the same four parameters in a slightly 
different form. The third transformation uses a vector to represent the 
rotation, and performs the transformation by means of purely vector and 
algebraic operations. 

These three types of transformation will be compared to determine 
which is most useful for machine computation. 
4.B Direction Cosines: Matrices and Tensors 


4.B.1 Direction Cosines 





The most familiar of the rotation transformations is the matrix 


transformation (Ref. 30; p. 95ff) 


Vy = CV, (4-1) 
where Vy is a column vector with components measured in an 


inertial coordinate frame. 
WE is the same vector with components measured in a 
body-fixed coordinate frame, and 
C is the square matrix of direction cosines of the inertial 
axes relative to the body axes. 

The first component of Eq. 4-l in vector notation is 

=i 


= (hy 


prep ty ba) po ea bs 
(4-2) 
where l,, is the unit vector along the inertial X axis, etc, Sinee the 


vectors involved are unit vectors, 


ln epiecee ae detalii (4-3) 
and in general, the elements of the matrix C are 
TT = cos (1; , Lap ere L; (4-4) 
The rate of change of the elements of C may be found by differentiating 
iq. (4-4) directly: 
ole heres re dale 
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jtl (4-5) 
or using a dot (*) to denote time differentiation, 


Ce aK 


1) )ataeal *j+2 C 


ij + 2 ~ j+1 (4-6) 
Note that ae is the jue component of angular velocity of the body -fixed 
coordinate frame measured in body coordinates. 
Eq. (4-6) may also be derived by differentiating Eq. (4-1) directly: 
V, = CV, + CV, » 
nae ° —] t . 
= + 
This is the Law of Coriolis in matrix form. If the term omen, be 


identified with w X V,, and the angular velocity matrix Q be defined by 


b 

0 W. “Wo 

oa ee 0 aT 
> Ria? (4-8) 

then 

SS , 2] 74 
sl. = 

w X Vy LV C CV,- (4-9) 


Thus, the differential equation for C is 


C = -CY, (4-10) 
which expands to give Eq. (4-6). 

There are nine direction cosines. Since a rotation involves only 
three degrees of freedom, there must be six constraints on these nine 
quantities. These constraints arise from the fact that during a rotation, 
the magnitude of a vector remains unchanged. This condition requires 
that 

Coe aeCs (4-11) 
where ae is the transpose of C 
and °° is the inverse of C. 


Another, more useful, form of the orthogonality condition is 
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ote = ce. =a (4-12) 


where I is the identity matrix. Thus it is possible to write the nine 
direction cosines as functions of any three of them, using six independent 
orthogonality relations from Eq. (4-12). Another solution involves the 
use of Euler angles. 
4.B.2 Euler Angles 

There are many ways of defining Euler angles; the definition 
shown in Figure 4-1 follows the most common usage (Ref. 30 p. 104ff., 
and Ref. 39, p. 23, for example). The first rotation is about the Z axis 
of magnitude An. The second rotation is about the displaced x axis, which 
becomes the line of nodes, by an amount Ay. The final rotation is about 
the z axis by an amount A. In example being used, the capital letter 
axes correspond to the inertially fixed frame, while the lower case axes 


are the body-fixed set. 





Fig. 4-1 Euler Angles 


In terms of these Euler angles, the direction cosine matrix is 
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cos Ay cos A, - gin A, cos Ay agin A, - cos A, gin A, - gin A, cos Ar cos A - ain Ay gin A, 
-sin A, cos A_ - sin A., cos A, cos A - sin A. sin A_ + cos A. cos A, cos A 
Zz z Z & z Z, z Z L Zz sin A; cos AL (4-13) 


gin A, ain AL "= cos A, sin Ay cos Ar 


The rates of change of the Euler angles in terms of the 
angular velocity of the body-fixed set measured in the body-fixed 
coordinates may be determined by resolution, and are (Ref. 39; 
p. 24) 

e ioe sin A, ee cos AL 
Z 1 sin Ay 2 sin Ay, 


>. 


2 cos A, ate 25 sin A, 


A, = ©, cot Ay sin A, - W5 cot Ay cos A, To, 
(4-14) 
This computation suffers from two disadvantages; trigono- 
metric functions must be determined, and a singularity exists 
at sin Ay = 0. The singularity corresponds to gimbal lock in 
a three-gimbal system (see Ref. 39, Figure 8-37, p. 25) and 
can be eliminated by including a fourth parameter in the 
description of the rotation, just as gimbal lock may be eliminated 
by the addition of a fourth gimbal. 
4.B.3 Tensors 
A third expression for the direction cosines is most 


conveniently cast in tensor form. It should be remarked here 
that since the transformation involves only Cartesian coordinates, 
no distinction between covariant, contravariant, and physical 
components of vectors and tensors need by made (Ref. 35, p. 170). 


The direction cosine tensor is 


, =cos A oi + (l-cos A) rivi- sin A Crh i (4-15) 


63 


where the rotation which carries the body-fixed axes into the inertial 
axeS is about an axis with direction cosines r; (the same with respect 
to either coordinate system) with magnitude A. Here ae is the 
Kronecker delta 
ag ea (416 
if if j 

Ciik is the cyclic tensor 

+1 if ijk positive cyclic (e.g. 123) 
=(-1 if ijk negative cyclic 

Oifi=j, j=k, ori=k (4-17) 


and the convention of summing on repeated indices is observed (Ref. 34, 


“ijk 


p. 96 for example). Using the tensor formulation, Eq. (4-1) is 


Vn = Sig%p3 (4-18) 


and Eq. (4-10) becomes 


Ci ik? kjm*m (4-19) 


4.C Four Parameter Representation 

The direction cosine computations are six-fold redundant; the Euler 
angle computations carry no redundancy but have a singularity. The use 
of four parameters to describe a rotation eliminates the singularity and 
diminishes the redundancy. There are two practically identical ways of 
describing a rotation with four parameters. These are Hamilton's 
quaternions (Refs. 31-33, 36, 38, 40) and Cayley-Klein parameters 
(Refs. 30, 40). 
4.C.1 Quaternions 

Hamilton developed quaternions as an extension of complex numbers. 
In representing a vector in two dimensions as a complex number, the 
unit real number 1 represents one direction, and the unit complex number 
i=N-l, the other. In analogy with this, Hamilton let the unit space 
vectors bei, j, k, with multiplication rules similiar to those for V-l1: 


%, 2 
oe: re ; = ne = -] ij = -ji =k, etc. (4-20) 


Thus, a space vector may be represented as V = iV, + jVo 4 kV., and 
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a general quaternion, with the unit real 1. the fourth direction, as 


q=atibt+jet+kd (4-21) 
Conjugation is defined by 

q =a-ib-je-kd (4-22) 
and length is given by 

0? = qq Sia +e + — a? (4-23) 


A simpler notation follows from representing a quaternion as a scalar 


and a vector: 
q = (x; ee Gs 1p) + jP5 7S Kp. (4-24) 
Thus; : | 
q = (A, -p) (4-25) 
and if q, = (A, py) and qo = 5; Po)s then 
SOAS er Bie SA yPoutepe ations) (4-26) 
Note that in general multiplication is not commutative.. Two special 


quaternions may be defined, the unit quaternion qa and the zero quaternion 
do’ 
Ske = (le 0) do = (0, 0) (4-27) 


Obviously 
99, ~> 42 = 4 494 ~ W945 ~ Uy (4-28) 
Further, if q 7 do? there exists a unique inverse or reciprocal 7 such 


-|_ -l 
thatqq =-q 94=q;: 


-| = l 3c 
q ~ 94 
lq! (4-29) 


Two relations follow from the definitions of conjugate and inverse: 


( ) * = KG _ ( i = “het 
449 So) hh 4149 BH (4-30) 


If the quaternions q, and q 1... defined by 
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ay = (2p Gee a 


(4-31) 
and related through the quaternion q having unit magnitude. 
q= (7) AP + p=] (4-32) 
by the transformation 
a, = 44,4 (4-33) 
then 
(a, V,) =[a,, Vit 2X6 Vi + 2p X.(0 X VI Py 


The scalars of q; and q, are equal, and the vectors are related by 


V, = Vj, + 20 XV, + 2p X (p X Vb). (4-35) 
Writing out Eq. (4-18) in vector notation, and using ras the unit vector 


along the axis of rotation gives 
V.=V, cos A+ (r- V,) r(l-cos A)+rXV_sinA 


1” ‘b b 
-V,+rx V, sin A +[(r-V,)r - (rt r)V,] (1-cos A) 
Seat rX V, sinAtr X (r X V,)(1-cos A) 


vie anon 4 cos : aie + 2 sin” = r X (r X V,,) 
(4-36) 
If X in Eq. (4-35) is identified with cos in Eq. (4-36), and p with 
r sin Zan then the transformation Eq. (4-32) corresponds to the rotation 
of V, about an axis r through an angle A. Writing q out gives 
q = (cos 4 r sin 4) - cos + sin4 (ir, + jro + kro) 
(4-37) 
Extending Euler's formula for exponentials of complex quantities 
Bb ft Ne the notation 
ia + JT 5 a3 kr) a cosa agg (ir, : irs ne krg) 


2 2 (4-38) 


with r try” ra” = 1, so that the transformation Eq. (4-33) may be 


written 
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A,. , AY, ’ 
; 5 (ir, oF jr 1 kra) ns (ir, + jr, + kr4) 
The quaternions representing successive rotations may be combined as 


follows: 
% oc 
Wo ~ Gp [ap 9K |) qa: (da dA ) 1 (da qa) 
2 l 1 2 Be | 2 (4-40) 

The areinQarara ene ae of finite rotations is evident from the non- 
commutativity of quaternion multiplication. 

The rates of change of the rotation quaternion parameters must now 
be related to the angular velocity of the body-fixed coordinate system 
with respect to inertial space. Let q be the quaternion representing the 
rotation that carries the body-fixed space into the inertial sett: at time 
t. The quaternion representing the opposite rotation is then q Let 
the quaternion representing the rotation carrying the body-fixed frame into 
inertial space at time|t + A t|be q'. The quaternion representing the 


rotation of the body-fixed axes during the time At is therefore q'q*: 


= [ (0, p) + (AA, Ap)] (A, -p) 


= (1, 0) 3 (Ad, Ap)(X,-p) 
(4-41) 


or 


AA Ake 2, : 
(cos aoe 1, sin = = (AA, Ap)(A, -p) (4-42) 


Dividing this equation by At and allowing At to approach zero gives 


AA 
ape cee 2 AA si A 3) 
At +0 | A A/2 “ADT At ~AA/2 Atl = (BEREAN, p) 
= lim | a : Z 


- AA A 


lim : : 
Since Ayn Yat at w, Eq. (4-43) becomes 
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1 — a a = - 
9 (0, w) = OS p A, -p) (4-44) 


or, post-multiplying by (),p) 


(1,8) =) (0,0) (A, 0) : 


whence (4-45) 


he -) Gp 
and 
p = 5 (e+ GXp) (4-46) 
4.C.2 Cayley-Klein Parameters 
The foregoing discussion provides a method for describing a 





rotation with four parameters. One is redundant; the redundancy is 
removed by the constraint that the quaternion q = (.,¢) have unit 
magnitude [ Eq. (4-32)]. This formulation suffers from the fact that 
quaternions are not at all widely known. To obviate this difficulty, 
Cayley-Klein parameters may be introduced; they retain the simplicity 
of quaternions while eliminating Hamilton's i, j, k. Defining the Cayley- 


Klein parameters by 


a =x - 1p, 

B = -Po - ips (4-47) 
ol TP 5 ~ 1P 2 

6 =X + 1p, 


and from them, forming the complex matrix Q, 


@ ae Bp 
Qu= vy 6 = CF s ok (4-48) 


gives a matrix which is unitary and has determinent +1: 


a 
* a -i 
Q =Q =Q 
RE tr ee 


ofc 3 9 (4-49) 
1Ql =aa +6B = SORE epaa io, SV ae ae ol 
a 
Here, Q@ is the complex conjugate of Q, eo is the transpose of Q,and 
Gia is the adjoint of @. When using Cayley-Klein parameters, vectors 


are represented as Hermition, or self-adjoint, matrices: 
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Vv = [| V3 Ne ae: (4-50) 
Vv, a iV. -V2 
Note that the trace of the matrix V is zero, and the determinant of V is 
the negative of the square of the magnitude of the vector V: 


ee ey, ae 
IVI = -(Vp+ V5 + V3) = -V re 


Using the matrix Q, a similarity transformation may be performed 
on the matrix V,, yielding another matrix Vy 
V, = QV,.Q" (4-52) 
I b 


Under similarity transformations, the Hermition property of a matrix is 
preserved, as are the trace and determinant. Therefore, Vy must be of 


the same form as V, [Eq. (4-50)], and further 


=> ; _ = 
-Ve=lv,l = lvl = -Vy 


Eq. (4-53) is the orthogonality condition. Thus the transformation 


(4-53) 


Eq. (4-52) is equivalent to the transformation Eq. (4-33). 

The differential equations for the Cayley-Klein parameters follow 
directly from the defining equation [| Eq. (4-47)] and the rates of change 
of the quaternion parameterslEq. (4-46)]. 

4.D Vector Representation 
Another useful representation of a rotation is the vector o 
o =r tan A (4-54) 
where r and A are defined above. In terms of the scalar and vector of 


the rotation quaternion 








ain nee 2 i 
cues aa oe aaa i ES 
cos= 2cos = 1+ cos = - (4-55) 
4. 4. 2 
The inverse relations are 
eee A 
Metros = 2 cos’ -l = g -l = 3 
2 4. 2A " 2A 
sec a sec Ve 
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7 Lo36 
; ae coe 
and : (4-56) 
p = (1 FX) Go eo oO 
Substituting these expressions for X and p into Eq. (4-35) gives 
—_—-— = l - o° 20 
V,=V,+ sedi: ee lous ecanr Boa) Ea) x Ti 
or _ 
wh= Vo + 4( - 0”) oxXV. + 8 oX (0X V,) (4-57) 
tn aiiaal edit attaaemmaliiiiall ; 


The vector representation of a rotation requires only three numbers; 
there is no redundancy. There is a singularity, however, at]Al = 27, 
corresponding to a full rotation. This singularity may be eliminated by 


defining another vector o': 


go' = _ 
oe 
7 (4-58) 
Gg) =——=—=(- fF) cot 4 (-r) tan Ci = 4) 
tan — 
4 
S (<a) tan (4-59) 


4. 

Since a rotation of A in one direction is equivalent to a rotation of 
2x ~-A in the opposite direction, o and o' describe the same rotation, and 
Eq. (4-57) holds with o' substituted for o. 

The differential equation relating the rate of change of 6 to the 
angular velocity of the body-fixed coordinate frame may be obtained by 
differentiation of Eq. (4-55): 


ee ee 
~ » oe 


Qi* 
i 
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_1,AS+6xXp , (3-6)5 
2 ea oN (1+ 





Ww +wX ot (Ww o)O] 
(4-60) 


w+twoXoat (wxo)X ot aaa) 








[ —s— & + & X OO! + (uXO') X O!] (4-61) 


4.E Comparison 


The preceding three sections of this chapter have presented several 
methods for transforming a vector from one Cartesian coordinate frame 
to another, and derived the differential equations for the transformation 
parameters in terms of the angular velocity of one frame with respect 
to another. These transformations and their differential equations are 
collected below for comparison. 


4.B Direction Cosines 


Transformation 

Vy = CV, (4-1) 
Direction Cosine Differential Equation 

C = -C2 (4-10) 
Six-fold Redundant, with Six Constraints 

Cle 27 (4-12) 


Euler Angle Differential Equation 


; sin AL cos A, 
A = eee ODOT 
VA, 1 sin Ay, Z sinA, 
A. =wW, cos A, oF Wo sin A, (4-14) 


L I. 


A = -w, cotA 
Z, 


1 L sin A, =B5 cot Ay, cos A, a We 
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No redundancy; Singularity at sin Ay = 0 


4.C Four Paremeters 


Quaternion Transformation 


Vy = Vi, + 2Ap X V+ 2p X (p X V,) (4-35) 


c Lie alanine ees : 
A= cos 5 p=rsin “} (4-37) 


Cayley-Klein Transformation 





Vi3 Vii- Ye A - Ip, “Pg - oy Vb3 Vor = Yba] J > * 3 Pg * 1, 52) 
Vy * iV yp “Vig Po- % Athy Voi * Vp2 “Vp | | -P2 * Ad - io, — 
Differential Equation 
a = a. p 
(4-46) 
o => (Xw +wXo ) 
Once redundant with one constraint 
Vins ae (4-32)(4-49) 
4.D Vectors 
Transformation 
V.=V_+4 1 - of oxXV_ + 8 -o X(o XV,) 
I b (1 + 02) b (1 + 92)? b 
(4-57) 
- - A 
o =r tan > (4-54) 
or 
72 Vo +4 ao XV, + g o' X (a' X V, 
I b (1 + ot2y2 b (1 + ot2y2 b) 
(4-57a) 
Fg Cias 
52 (4-58) 


ee 


Differential Equation 


2 aie aun 

G = orl 5 w+wXot+ (wXo) X o] (4-60) 
a ete ? 

Gt => [SSS + OG! + (GKE")X 0") (4-61) 


These transformations will now be examined with a view to determining 
which is most suitable for machine computation. The basic criterion to 
be used in this determination is the size of computer required. Further, 
since the information being dealt with is incremental, incremental 
computation will be the basis of decision unless general propose computation 
seems particularly appropriate. 

The most straightforward computation employs the direct direction 
cosine formulation. The transformation requires nine incremental 
integrators, since each component of the vector to be transformed has 
magnitude AV. The generation of the direction cosines requires eighteen 
incremental integrators. If Euler angles are generated, then the sine 
and cosine of each angle is also required; further, the sines and cosines 
must be combined to give the direction cosines as well as the factors in 
the differential equations. The use of incremental devices requires 
forty-six integrators. This does not handle any problems raised by 
the singularity at sin Ay, = 0 

Twelve incremental integrators are necessary to generate the four 
parameters used in the quaternion and Cayley-Klein parameter representation. 
If the direction cosines are formed from these parameters, an. additional 
twelve integrators are required. Nine more complete the transformation. 
Thus, thirty-three incremental integrators are required in all, as opposed 
to twenty-seven for the direction cosine computation. If the formal 
Cayley-Klein similarity transformation is used, the first multiplication 
can be performed by twelve incremental integrators. An arithemetic 
computation seems best for the second matrix multiplication; it requires 
twelve multiplications and nine additions. 


Twenty-one incremental integrators are required to generate the 
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three components of the vector used to represent a rotation; the vector 
transformation requires twenty-five more. In addition, an arithmetic 
section is required to perform the branching operation used to remove 
the singularity. 

The foregoing discussion indicates that the direction cosine 
computation using a total of twenty-seven integrators is most compact. 
For this reason, this formulation will be the only one considered in 
designing the transformation computer. 
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CHAPTER V 


TRANSFORMATION COMPUTER 





b.A. Introduction 


In Chapter IV, the various ways of representing coordinate ro- 
tations were derived and compared to determine their suitability for 
machine computation. As aresult of this comparison, the matrix trans- 
formation was chosen to be instrumented, The rotation matrix is to be 
generated by direct computation of the direction cosines, The matrix 


differential equation to be instrumented is 
C= -C2 (5-1) 


where C is the matrix of direction cosines, a dot (-) Over a quantity 


denotes time differentiation, and 92 is the angular velocity matrix, 


Q = “We 0 W; (5-2) 


The angular velocity is resolved into body axes. 
Eq. (5-1) represents nine separate equations; there are three | 
sets of the form 
Csi = age 
oe eb iecales el 


Cig = Ci - Cio) 


Multiplication of Eq..{5-3) by dt and integration of the result yields 
Ci, = J Cj9d83 - J Ci3d0, 
Cig= J C40, = J C208, 


5.B Algorithms Available 


Equations (5-4) show that each direction cosine is the difference 
of two integrals, each the integral of a direction cosine with respect 
to rotation about a body axis. Since the gyro outputs are equal incre- 
ments of rotation about body axes, the integrals may be approximated 


by quadrature rules using equal A® increments. 
5.B.1 Quadrature Rules 


Almost any standard text on numerical analysis (Refs. 41, 
46-48, 52, 54, for example) gives the more common quadrature 
rules. Since only the integrand will be available, no rules employing 
derivatives of the integrand, such as Taylor's series, will be con- 
sidered. The most common quadrature rules using present and past 


values of the integrand are the rectangular rule, 


@6(n) 
C(@)d8 = AGC[O(n-1)] 
" Q(n-1) (5-5) 
the trapezoidal rule, 
O(n) 
(crayae = 88 (CfO(n-1)]+C[@(n)] } (5-6) 
Q@n-1) 


and Simpson's one-third rule, 
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Q(n) 
c(e)de = 2° (cle(n-2)|+ 4C[(n-1)]+C[9(n)] } (5-7) 
(n-2) 


The first is an open rule; that is, the value of the integrand at 
the present time does not appear. The remaining two are closed 
rules, Each integration rule fits a (p-l) degree polynomial through 
the p points of the integrand appearing in the formula; for this reason, 
these rules are called (p-l) order polynomial rules. The rectangular 
rule is called a zero-order rule; the trapezoidal rule and Simpson's 
one-third rule are first- and second-order rules respectively. 

At each step, the approximation to the integral is in error by 
some amount. This error, arising from not considering higher 
derivatives or differences, is called truncation error, For the three 


rules above, the truncation error made at each step is (Ref. 47, p. 73). 


Rectangular EF. - 0(A0%) 
Trapezoidal E, = 0(Ae?) 
Parabolic  E,, = 0(A0°) 


This notation means that for the rectangular rule, the truncation error 
at each step is on the order of Ae7, The term "on the order of'! 


means that the limit of E,/Ae" remains finite as Ae? approaches zero. 
9.B.2 Influence of Size of Step 


The rather fast pulse rate required in the gyro loops suggests 
that using higher order integration rules with larger increments 
would reduce the speed requirement on the computer as well as de- 
creasing the truncation error. For example, if rectangular integra- 
tion were used with a step length AQ, the truncation error at each step 
would be 0(Ae7), If trapezoidal integration were used with a step 
length of 5 A®, the error at each step would be 0(125 A?) which is 
Significantly less than the rectangular error, This increase in step 
size permits reduction of computer speed by about a factor of two. 


However, there is an irreducible uncertainty in the gyro information 


77 


of 0(A07) in each pulse, so that the apparent increase in accuracy may 


not, infact, be attained. As is shown in Chapter III, Section 3.B, 
nT, 

each output pulse of a gyro is {usa , the integral of one component 
ov) Le 

of angular velocity resolved into body coordinates. Since addition of 


rotations is not commutative, 


nT, nT. 

- 2 — s 

\ dt -#i (aya +j \wodt+k . \wedt (5-8) 
ola Dy He ba Goa se eet Iam (n-l)T. 


in general. In one sampling period, the angular velocity could be 
equal about all three axes, or there could be successive rotations 
about all axes in any order; each input gives the same output. 
Eq. (5-8) adds these rotations as though they were vectors. The 
error thus made is 0(A07) (Ref. 34, p. 97). The conclusion to be 
drawn, therefore, is that rectangular integration of Eq. (5-1) is the 
most suitable algorithm because the angular information is available 
in rotational increments. Note that this comment applies to any 
digital computation where the rotations are handled in this way; it 
is not restricted to delta-modulated gyro information. In special 
cases, however, the use of higher order integration rules to reduce 
truncation error may be of some advantage. 

The equations to be instrumented to generate the direction 


cosines are three sets of the form 


AC, = Cj9A8, - C;,A85 
DGioe= Cg ie CAO (5-9) 
ACi3 = Cj A8_ - CjoA8 


9.C Digital Differential Analyzers 


The computer chosen to generate and perform the coordinate 
transformation must be fast and accurate. The accuracy require- 


ment, with present equipment, eliminates analog computation from 
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ee See 


(a) 


Register Arrangement 


AX 


EDK 


AZ 


AY Y 


(b) 


Schematic 


Fig. 5-1 DDA Integrators 


AY 


consideration; the speed requirement, based on the necessity for 
processing each pulse as it comes from the gyro, calls for a com- 
puter whose cycle time is no greater than the sampling period of the 
gyro loop. The sampling period, as determined in Chapter IU, 

Sec. 3.B.4, is the quotient of maximum angular velocity and angular 
increment. 

A very fast and compact type of digital computer, and one 
particularly well suited to the solution of Eq. (5-9),is the digital 
differential analyzer (DDA) (Refs. 43-45, 53). 

The DDA has been described as an analog computer in digital 
form; it has the advantage found in the Bush differential analyzer 
(Ref. 42) of not requiring time to be the independent variable of 
integration. The basic element of DDA computers is the DDA 
integrator, usually shown diagrammatically as in Fig. 5-l{a). The 
blocks labeled R and Y are storage registers. The line labelled 
AY carries pulses representing a change (positive or negative 
according to the sign associated with the pulse) in the least significant 
bitin Y. The AX line carries pulses representing equal increments 
in X; when a AX pulse occurs, the contents of the Y register are 
added into the R register or subtracted from the R register, depending 
on the sign associated with the AX pulse. The contents of the Y reg- 
ister are left unchanged by this operation. The AZ line carries over- 
flow pulses from the R register representing a change in the least 
Significant bitin Z. The AZ pulses may be used as AX or AY pulses 
in another integrator. The name "integrator" is applied to this unit 


because it computes 


AZ = YAX : (5-10) 


This is an increment of rectangular integration of Y with respect to 
X [See Eq. (5-5)]. 

Each of the variables X, Y, and Z represents a physical 
quantity to some scale, If the Y register is n bits long, so is the R 
register. The least significant bit in the R register, AR, is related 
to AY and AX by 
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AG; 


ACia| Ci2 





Fig. 5-2 Arrangement for DDA Solution of Equation (5-9) 
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AR = AXAY (5-11) 


The overflow logic for the R register causes an overflow pulse when 
the content of the R register exceeds AZ. The simplest kind of 


overflow logic scales AZ so that 
AZ = 2°AR (5-12) 


Thus, the AZ overflow occurs when the R register fills up. AZ is 
then related to AX and AY by 


AZ = 2UAXAY (5-13) 


The usual schematic representation for a DDA integrator is 
shown in Fig. 5-1 (b). 

The solution of Eq. (5-9) requires six incremental integrators 
for each value of i arranged as in Fig. 5-2. Note that no order of 
performing the additions is specified. 

There is feedback in this arrangement; hence there are scaling 
consistency equations to be satisfied, just as in the case of analog 
computation. The integrators are identica]; thus consideration of the 
loop with the AB. integrators leads to 


nm 2 Zn _ : 
A8,AC.,2 = A8. AC.52 - AC.» (5-14) 


or 


A8, = an (5-15) 


by repeated application of Eq. (5-13). A similar result holds for 

A8, and 485. Thus for simplest overflow logic, A@ should be chosen 
as a negative integral power of 2. The value which A®@ must not exceed 
may be found from system specifications, andn, the length of the Y 
and R registers, is determined. Requiring that the scaling of the 
velocity increments be the same in both coordinate frames yields 

INC 2 


82 


o.D Analysis 


With the form and function of the transformation computer de- 
termined, an analysis of its operation must be made. Three points 
are usually considered in the analysis of a computer: dominant mode, 
truncation error, and round-off or quantization error. These will be 
taken up in turn after a frequency domain model for the computer has 


been developed. 
o.D.1 Frequency Domain Model 


Computers operating in real time have been characterized as 
sampled-data systems and transfer functions derived for them 
(Ref. 51). Since the equations that the computer is solving are dif- 
ference equations, they may be transformed in a way similar to that 
used in transforming time sequences or samples of continuous func- 
tions (Ref. 61, Ch. 1). 

The equations describing the computer operation with 


W) = Wo = 9, We # 0 are 


C,,(n) = C,,(n-1) + A8,C..(n-1) 
(5-16) 
C,(n) = C,9(n-1) ~ A®.C.)(n) 


_. A 


Let Z. =e Ws, where s is the Laplace complex frequency 


variable. Then, if 
_ n 4 
CA) = 2M (5-17) 


Eq. (5-16) may be transformed into 





A9,Z. 
C5,(23) =————_ C9 {z3) 
1 - Ze 
(5-18) 
AQ 
Ci9(Z3) = - C5 ,(23) 
l-Z. 
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if W) = 2Wo, We = 0, the transformed set of equations is (see App. C, 


Sec. C. 1.5) 


2 2 
Ae» 1 +=. 


aan eae’? 
Cia a: SS Ci3(2) 
("27+ 2) oe ee 


A8z 
= C;,(z) | (5-19) 


Cj 9(z) = 











mN: Ae 
Ci3(2) = 5 Cy 2 - Cio(z) 
l-2z l-z 


rr) 
Mo > 


This set of equations may be well approximated by 


beeen -s 2 
where z =e a and z =e nea eS 


A8z 











C, (2) ai C; (2) 
l-z 
CF FS ez) (5-20) 
l-zZ 
mn: AQ 
Ci3(Z) = Cj ,(z) - —— C;9(z) 
l b-z 
Z 
Ae aX 


If the variables Zy and Zo are defined by Z) = ere W) and Zo =e Wo s 


then Eq. (5-20) may be written 








7 A8z 
C,,(z) ae C;3(z) 
] - Zo 
— Aoz 
C; (2) - C,3(z) (5-21) 
] - Zy 








Cig(z) = 9° c,2) - 92 cz) 
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where z represents the delay between the cycling of integrators. If the 


components of angular velocity are related by 


kw, =fWo = MWe = w (5-22) 
SS A®8 
where k, @, and m are integers, and if Zz, =e a4 ,» then the Zs! s are 
related by ya 
i - a = Ze = 7 kim (5-23) 


Here again z represents the delay between the cycling of the integrators. 
With the above relations between angular velocity components satisfied 


the set of equations [cf. Eq. (5-9)] 
C,,(n) = C,,(n-1) + A8,C,5(n-1) - A85C,; 3(n-=1) 
C,.(n) C,9(n-1) + A®, C;3(n-)) - A83C.,(n) (5-24) 


C,(n) C,4(n-1) + A0,C,(n) - A0,C;o(n) 


may be approximately transformed into 














_ Avz A8z 
C.,(z) “a C. (2) se C; (2) 
=e aS 
_ Adz AQ 
- -Z 
3 
_ Ae A8 
ae, aml 


With Ww) = 2Wo, We = 0, this set is identical with Eq. (5-21), which 
is a good approximation to Eq. (5-19). Similarly with w, = Wo = 0, 
We #0, Eq. (5-21) becomes Eq. (5-18). Note that if Ws 0, the terms 


in Zs are taken as zero. 
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Eq. (5-25) may be represented by a signal flow graph as shown 
in Fig. 5-3. When the angular velocity of the system is specified, the 
substitutions indicated by Eq. (5-22) and Eq. (5-23) may be made, 
For example, with W) = 2Wo = 8W9, k =1, @ = 2, m = 3, and the appro- 
priate signal flow graph is shown in Fig. 5-4. 

The use of Mason's method for reducing signal flow graphs 
(Ref, 60, Ch. 4) shows that the characteristic equation for the com- 
puter under these conditions is 

2 2 2 


5 
14 AR 2% , A492 4 492 z AQ z = (0 (5-26) 


2 9 
(l-2)* (-2*) (-29)° (-2%)a-2%) 


and, in general, the reduction of the signal flow graph in Fig. 5-3 


leads to the characteristic equation for the computer 


2 2 2 3 
1 +49 2 4 AVzZ + 49 2 _ As z(l - z) 7 
l 2 3 
(5-27) 


With the node-to-node transfer functions determined, the remainder of 


the analysis can proceed. 
5.D.2 Determination of Dominant Mode 


The various normal modes of the transformation computer are 
determined from the roots of the characteristic equation, Each root ¢ 
leads to a time sequence of the form aaa with coefficients chosen to fit 
initial conditions, If € is complex, its conjugate will also be a root, 
and an oscillatory term will be present. If |I€| >1 then the mode will 
decay, while if |¢| <1 there will be exponential growth. 

The initial conditions as the varicus nodes appear as step inputs; 
for this reason, the roots near z = 1 will provide the dominant modes, 
while other modes, which occur near the various no roots of unity, 
will have magnitude A989. Thus the dominant mode of the computer 


comes from the roots near z = +l. 
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Fig. 5-3 Signal Flow Graph 


Cil 





Fig. 5-3 Signal Flow Graph for Computer Fig. 5-4 Signal Flow Graph for 


WW, = QW. = 3wW 


l 2 3 





Fig. 5-5 Analog Computer Flow Graph 
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The magnitude of the roots and their associated modes may be 
found by root-locus techniques. Since the loop gain, Ne. is so 
small, linear approximations must be used in conjunction with the 
graphical technique... 

In Appendix C, the analysis of two situations is carried out 
using root-locus techniques. For the case of angular velocity ratio 
two to one, there are extraneous roots near the square roots of unity 
not +l. In the case of equal rates on all axes, there are only roots 


near +l. 
5.D.3 Truncation Error 


Truncation error is the error made by not considering higher 
derivatives or differences when computing the solution of a differ- 
ential equation. Itis the difference between the exact solution and 
the computed solution using computer registers with infinite length. 

Assumption of the same angular velocity relationships as above 
[Eq. (5-22)] with the added restriction that the angular velocity be 
constant, permits the solution of Eq. (5-1). With Eq. (5-22) satisfied, 


Eq. (5-1) becomes three sets of the form 


~  . W W 

Ci tem ~9i3 7 

S Wo ce, W 8 
Cio = Cig - Cam (Sag) 
ry ‘ ag P w 

Cig = Cig - Cia k 


Applying a Laplace transformation to Eq. (5-28) yields 


C,,(s) == is Cio(s) - = Cy 3(s) + Cy) (0) 
s [m g 


Ciols) = 


wm |e 


W W 
F C,,(s) - = C,j(s) + C910) | (5-29) 


_lfw W 
C,,(s) = F C,,(s) a ‘ C,(s) a C53 (0) 
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where the C,,80) are initial conditions. Converting this set of equations 
to a signal flow graph leads to Fig. 5-5. The initial conditions are taken 
as inputs at appropriate nodes. Flow graph reduction may be used to 


derive expressions for the signals at each node, For example, 








2 Co(0): 9 CEO) C.4(0) Co) 
C 169) ? +a] ee saa | stu) 8 + i2 
Saale) 2 k g m km ke 
(k? +07 + TCE See ne Aer | 
(k- 2pl> See 
(5-30) 


These expressions can be transformed into the time domain and the 
exact solution to the differential equations found. 

The solution to the difference equations instrumented by the com- 
puter may be found by taking the z-transform of the solution, as deter- 
mined by root-locus methods, and writing down the corresponding 
time sequences. There will be truncation error from three causes; 
the roots near unity may not have the proper argument, they may not 
be on the unit circle, and there may be other roots. The error in 
argument of the root near unity makes the oscillations of the direction 
cosines the wrong frequency. If these roots are not quite on the unit 
circle, the oscillations will not be of constant amplitude as they 
should. The other: roots may lead to exponential or oscillatory 
terms; in any case the size of the corresponding mode may be deter- 
mined. 

When finding the z-transform of the computer output, the last 
root-locus construction, to find the roots of the computer character- 
istic equation, deals with a loop that has a gain of Mee Since Ae? 
is very small, the construction cannot be carried out exactly, but 
linear approximations may be made analytically. These approxima- 
tions are sufficient to make the truncation errors predicted inexact. 
However, the order of the errors predicted is correct and the growth 


rate is correct, but an exact error analysis is not provided. A more 
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Fig. 5-7 Quantization with Half-Grain Hysteresis 


90 


nearly exact error analysis may be had by finding the roots of the 
characteristic equation exactly, and finding the magnitudes of the 
modes by root-cocus techniques, but the present analysis provides 


enough information to form a basis for engineering decisions. 
5.D.4 Quantization Error 


Quantization, or round-off as it is sometimes termed by 
numerical analysts, occurs whenever a number is assigned toa 
physical variable. Through the quantization process the value of 
the variable is represented as the nearest integral number of the 
basic units in which the variable is being measured. If this basic 
unit be q, the quantization process may be represented by the input- 
output relationship shown in Fig. 5-6, 

While this is the basic quantization process, the overflow logic 
of many computers causes hysteresis in the quantization process. 
This hysteretic quantization has the characteristic shown in Fig, 
o-7(a). The hysteresis may be represented by feedback around the 
basic quantizer as shown in Fig, 5-7(b). Note that the quantizer has 
grain size 2q. The z in the feedback loop represents a delay of one 
Cvele tinmesoisthe Computer. 

Since quantization distorts the signal, although in a deter- 
ministic way, the distortion may be thought of as noise, Widrow 
(Ref. 55) shows that if the amplitude of the input to the quantizer 
ranges over more than three quantization boxes, the quantization 
noise, to a good approximation, may be treated as an additive first 
order random process, q(t). The probability density distribution 
of the quantization noise is shown in Fig, 5-8. The expressions 
for the autocorrelation function and the power spectral density of 


quantization noise are 


2 
(rat) 
fad 15 


(5-31) 
2 


ra) (z) = a 
qq 12 
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Fig. 5-9 Representation of a Quantizer by Independent Additive Noise 
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Fig. 5-9 shows the representation of a quantizer as independent additive 
noise. 

Although one attempt was made‘to apply Widrow's theory to 
DDA's in real time (Refs. 49, 50), some difficulty was encountered 
in applying it to interconnected DDA's., 

The appropriate signal-flow graph for the computer with non- 
hysteretic quantization is shown in Fig, 5-10, and Fig. 5-1ll shows the 
quantizers replaced by their equivalent gains and additive noise. 

If hysteretic quantization is used, a self-loop with gain +z/2 
must be included at each node, With integrators cycling at different 
rates, the z used in the self-loop must correspond to the fastest basic 
rate: the z of Eq. (5-23) must be used. 

Care must be taken to insure that a zero value of the variable 
actually corresponds to the zero of the computer. In many cases a 
bias of 54 is introduced because of an ambiguity of initial condi- 
tions (Ref. 49, p. 51). 

Using the statistical model of quantization and the sampled- 
data model of the transformation computer, an expression for the 
power spectral density of quantization noise may be developed. For 
example, with W) = ZW = SW, the power spectral density of quan- 


fizacion in Ci is 


2 
2 2 9 
AOA , AOy ye Oe , S82 ,__Ad*2 
® (z) = ee T= ws? (= as 2p (= 2) Ke- 
il AG 
2 2 y) 3 12 
ye SO 2 ,AQ ,AGZ - Az 
WPS ee?) ea sh 2") 
(5-32) 


with similar expressions for power spectral densities at other nodes. 
Note that for the case of the computer, the quantization grain size 
is AQ, 

The mean square quantization error may be obtained by 


forming the sum of the squares of the sequence represented by the 
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Fig. 5-ll Computer with Quantization Noise 
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transfer function in Eq. (5-32) and multiplying this sum by the mean 
square of the quantization noise (Ref. 61, Ch. 7). Since in the case 
of this computer the roots of the characteristic equation have magni- 
tude close to unity, this mean square value may become very large. 
Thus, it is appropriate to consider only a finite number of terms 

in the sum-of-squares series corresponding to a finite number of 
steps in the computation. With the series truncated after a finite 
number of terms, an estimate of the error may be obtained. An 
example is shown in Appendix C, Section C.3. It is found that 
quantization error is much smaller than the truncation error; the 
number of steps required to make the rms quantization error equal 
to A@ is 0(A8~%). 

In the same section (Section C. 3), an expression for 
hysteretic quantization error is derived. It is found that the mean 
square value of hysteretic quantization error is greater than that 
of non-hysteretic quantization by a factor of four because of the 
larger quantization box, For this reason, non-hysteretic overflow 


logic should be used in the transformation computer. 
9.D.5 Reversibility 


Up to this point, the analysis has assumed that all the A8@ 
pulses have been of the same sign. Of course, this will never be 
the case, and some examination must be made of the effects of 
negative pulses. 

In the sine-cosine generators shown in Fig. 5-l2, first one 
variable is updated, and then the other variable is updated. Thus, 


Fig. 5-l2(a) implies the difference equations 


cos (nA8®) = cos [(n-1)A8] - A® sin [(n-1)A9] 
(5-33) 
sin (nAQ) = sin [(n-1)A98] + A9 cos (nA8) 


If a positive pulse is followed by a negative pulse, there should result 


cos [(n+l)A®] = cos [(n-1)A®] (5-34) 
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Fig. 9-12 Reversibility of Sine-Cosine Generation 


The use of Eq. (5-33) for processing both pulses yields 
cos [(n41)A8] = cosuffn-1)A0] + AQ“ cos [(n-1)A0] (5-35) 


If positive pulses are processed according to Fig. 5-12(a), however, 
and negative pulses according to Fig, 5-12(b), Eq. (5-34) is satisfied, 
the computation is completely reversible, and the truncation and 
quantization errors for the two steps cancel each other, 

This idea may be extended to the direction cosine computation, 
In Fig. 5-3, all the pulses are assumed positive, and the order of 
processing the pulses updates CH first, then Cio: then C53. Since 
the computation is made up of a number of sine-cosine generators 
interconnected, the process can be reversed by moving the z in the 
numerator of each loop transfer function to the dther branch; this 
may be accomplished by updating the variables in reverse order 
for negative pulses, 

There are three pulses, each with two possible signs, which 


drive the computer. This gives 2° = 8 possible combinations. 
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There are only three direction cosines ina set; the number of possible 
orders of updating them is 3! = 6. Thus one pair of combinations 
cannot be reversed. Using the updating sequence 1, 2,3, for all posi- 
tive pulses, the computer-~cannotsreverse . A®, ahs A8,-,A8.-; 
and A8.-; A8,t, AG st. With another sequence representing all pulses 
positive, another pair could not be reversed; however,any combina- 
tion of two, as well as any single pulse, can be completely reversed. 

The non-reversibility of this computation causes a drift of 
0(A0%) for each step reversed. Unless the computation is operating 
with marginal accuracy, the added equipment required to change the 
order of updating the variables is probably too expensive, in terms 
of weight, power, and cost, to pay for the additional precision 


achieved. 
o.E Processing AV Pulses 


At each sampling instant in the accelerometer loop a AV 
vector, resolved into body-fixed components,appears at the accel- 
erometer outputs. If a two-level relay is used in the loops, there 
will be pulses at each loop output at every sampling instant; ifa 
three-level relay is used, output pulses will occur only at certain 
times. The AV vector may be resolved into inertial coordinates 


by processing it with the direction cosine matrix: 


AV, = CAV. 


I i. (5-36) 


Since the magnitudes of the AV pulses are the same, Eq. (5-36) may 
be realized by the arrangement shown in Fig. 5-13. The AV, pulses 
will be accumulated in Vy registers. Errors also arise in this 
computation, and they may be determined by the methods outlined 


above. 
o.F Summary 


In this chapter, the design and analysis of the transformation 
computer was completed. Various quadrature rules were examined, 


and the rectangular rule was chosen as the most appropriate for 
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integrating the direction cosines [Eq. (5-9)]. With the angular 
information available in A® pulses, digital differential analyzers 
seem the most suitable form of computer for the computation 

(Fig. 5-2). A frequency domain model for the computer was de- 
veloped (Fig. 5-3) and analyzed. The truncation error was estimated 
by developing an exact solution for the differential equation, and 
comparing it with the computer solution as indicated by root-locus 
plots. Quantization was treated statistically and it was found that 
hysteretic quantization caused a greater error than non-hysteretic. 
The integrators were shown to drift by O(ae?) for each computation 
step reversed. Finally, a method for transforming the AV vector 


was established. 
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CHAPTER VI 


SYSTEM DESIGN 


6.A General Considerations 

The previous chapters have described and analyzed the component 
parts of the system; the design parameters must now be chosen so 
that the system operates in an optimum manner to carry out the 
mission assigned. ‘T’'wo decisions may be made without reference to 
the mission which the system will perform: whether the instruments 
should limit cycle, and how the clocks in the gyro and accelerometer 
loop should be related to each other. 

6.A.1 Limit-Cycling Instruments 

In Chapter V, Section 402.5, it is shown that for each step reversed 
in the direction-cosine computation, the integrators drift by an amount 
0(A0*). For this reason, it is desirable not to have the gyros limit 
cycling. With a three-level relay in the gyro loop, the torquer current 
may be turned off during the quiescent periods, thus saving a great 
deal of power. Turning the current on and off, however, makes 
stabilization of the current more uncertain. 

Since the velocity pulses are being processed by DDA's, there 
may be a drift of 0(A9AV) in accumulated velocity for each reversal. 
Thus, the arguments above indicate the desirability of eliminating 
the limit cycles in the accelerometer loops. 

6.A.2 Instrument Clock Relationships 

In general, the clock rate of the accelerometers will be slower 
than that of the gyros; thus the direction cosine may change many 
times between A V pulses. Processing of the A V pulses must 


occur between changes of direction cosines; this requirement may 


be satisfied either by providing buffer storage for the A V pulses, or 

by choosing the clock frequencies and phasing of the gyro and accelerometer 
loops so that the A V pulses will always occur at the same point in the 
computer cycle. The latter choice is more attractive, since it 

eliminates the need for buffer storage and the consequent delay in 
processing information. 


6.B Basic Design Parameters ; 





The quantities needed to design a guidance system are the inputs, 
or maneuvers it will experience, the allowable error at the target, and 
the time of flight. The input specifications must include the maximum 
acceleration, velocity, and angular velocity which the guidance system 
is expected to undergo. From these data, the basic instrument and 
computer parameters may be determined. 

With increment size and clock frequency specified, Chapter III and 
Chapter V indicate methods for determining the remaining system 
parameters. In the two sections that follow, computations are carried 
out to determine clock frequency and increment size for typical aircraft 
and missile missions. Numbers representing current performance are 
substituted into the expressions for increment size and clock frequency 
to provide a notion of the order of the magnitudes involved. 

6.C Airplane Flights 

~~ Airplane guidance systems are position control systems; their 

task is to guide the aircraft to a specified position. Thus, errors in 
position at the target are errors in the controlled quantity. To determine 
specifications for the instrument increment sizes, it 1s appropriate to 
analyze a ‘worst possible'' case, and equate the error made to the 

error specification. 

An airplane is essentially a constant speed vehicle; thus, the worst 
possible case for the generation of error is a flight at maximum velocity 
with the direction of the velocity in error by > Ad, and the magnitude of 
the velocity in error by > AV. If the allowable error in position a is 
the same for gyro and accelerometer quantization, the appropriate 


calculations are 
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al er: 
Pe BUN wan Ts 
and 
l EE. 
AV = 
a r (6-1) 


where Ty is the time of flight. Thus the required increment sizes are 


AG = + omamppensdias is 
V T 
Hvac et 
and oR (6 -2) 
AV = TT 


Note that A@ should be chosen as a power of 2. The clock frequencies 


are 
W max 
tous AG 
and (6 -3) 
—_ a max 
a AV 


Specifications which represent a typical high performance airplane are 


7 7 Z 

i el a. ft/sec 
Ww = 3 rad/sec 

max 

(6 -4) 

V = 600 kt 

max 
E/T, = .1kt 


Substitution gives 
AG = 333 urad 
AV= .2kt = .333 fps. 


Since A@ should be chosen as a negative integral power of two, take 


ag = 27)4 = 244 urad 
Then the clock frequencies are 

7 =12.3 ke 

f. = 482 cps 


Since it is desirable to have f. a subharmonic of te set, 


. 12300 _ 
4 : = 492 cps 
The complete set of design parameters is 
AO = 244 urad. 
f —alaea ke 
6 (3-5) 
AV = .333 fps 
f = 492 cps 


The period of the gyro loop is 61.3 microseconds. With 18 additions 
to be made in the transformation computer in this time, an add time of 
4.5 microseconds per word is required. Since a first-rate arithmetic 
computer with central control would be required in toto to accomplish : 
this task, a DDA is clearly the best machine for this job. It will be 
smaller, faster, and easier to build. 

For serial addition of a word with thirteen bits (twelve plus one for 
sign) in four and a half microseconds, a bit rate of about three megacycles 
is required. This performance is available in present DDA's. 

With increment sizes and clock frequencies fixed, the remaining 
parameters may be quickly determined. 

6.D Missile Flights 

In contrast to an airplane guidance system, a missile guidance 
system is essentially a velocity control system, and a specified error 
in point of impact must be translated into an error in velocity at thrust 
termination. Relations between increment sizes and velocity error at 
cutoff are determined and translated down range using linear perturbation 
techniques to find impact error. 

The relationship between trajectory parameters and point of impact 
is called the hit equation (Ref. 65). The perwrbation factors are called 
range derivatives and are obtained by differentiating the hit equation. 
The velocity range derivative, 56R/6V, is the one of interest in the 
following computation. 


The worst possible case to be analyzed here is where the velocity 
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is in error by - AV in the direction of maximum range oe ve: and 
where the error due to gyro quantization is a velocity error of =- V5 49 
in the direction of maximum range derivative. The origin of the first 
error ig rather obvious; the origin of the second is not. 

Suppose that as each “a pulse is processed, the direction cosines 


aren error by $ AC -= 349, then the error at each step is > ABAY. 


Summing this over the total velocity accumulation yields 


E(v) = = 4 agav = tagzav = tadv 
2 2 co 


2 
(6-6) 
The appropriate calculations for increment size are 
l oR _ 
p OU, OV ae 
and (6-7) 
1 oR _ 
5 SV oy a 
or 
25 
Se SAR 
CORT, | 
and 
(6-8) 
IN 
AV = P 
oR 
VU 
The clock frequencies are then 
W 
p = —max 
aan (6 -9) 
max 
a AV 
Typical specifications for an 1800 nautical mile missile are: 
Z 2 2 
Anax 7 lg=322 ft/sec 
WwW = .2 rad/sec 


max 
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V__ = 20,000fps 
CO 


FE = 300 ft (6-10) 


Substitution gives . 
A9 = 33.3 uradg 


AV = .333 fps 
Since A@ should be a negative integral power of two, take 
NG = 2iae 


= 30.5 urad; 
l 


The clock frequencies are then 
f = 6.55 Ke 
g 


and 
fs = 484 cps 
Since the accelerometer clock frequency should be a subharmonic 


of the gyro clock frequency, take 


_ 6950 _ 
f. 7 es 904 cps 


The complete set of parameters is 
Ga 30,50 rad 





fee = 6.00 ke 

g 
AV = .333 fps (6 -11) 
f = 504 cps 


a 
As in the case of the airplane mission, the computer must be 


quite fast, requiring an add time of about eight and a half microseconds 
per word. Again, the logical choice for the computer, is the DDA, 
since an airborne computer must be compact. Further, the advantage 
of flexibility associated with an arithmetic machine would not be 
realized because of the fast add time required. For a serial DDA to 
perform eighteen additions of words sixteen bits long in the computer 
cycle time, a bit rate of about two megacycles is required. This 


performance is available in present DDA's. 
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CHAPTER VII 


SUMMARY AND CONCLUSIONS; SUGGESTIONS FOR 


a me a ee 


FURTHER STUDY AND DEVELOPMENT 


7.A Summary and Conclusions 


> ——— 


This work has been concerned with a broad theoretical analysis of 
gimballess inertial reference equipment using delta-modulated instruments. 
The equipment considered consists of three body-mounted, delta-modulated, 
Single-~degree-of-freedom integrating gyros, three body-mounted, delta- 
modulated, single-axis pendulum accelerometers, and a transformation 
computer. The delta-modulated gyros deliver output pulses representing 
increments of rotation and the delta-modulated accelerometers deliver 
output pulses representing increments of veiocity This representation leids 
naturally to the choice of an incremental computer io perform the coordinate 
transformation. An examination of performance required from the equip- 
ment concludes the investigation. Following is a brief summary of the 


Specific points covered, and the specific conclusions reached. 


1. Gimballess inertial reference equipment is theoretically feasible, 
and the required computation can le done with DDA's currently available. 
The major obstacle to the realization of practical gimballess inertial 
reference equipment is the lack of inertial instruments specifically 


designed for this application. 


2. Delta modulation is a pu!se modulition system in which each out- 
put pulse represents a fixed increment in the input, or one of its integrals 
or derivatives. The analysis of a deita-modulated loop proceeds ina 
straightforward manner, using a moding plot to determine the modes of 
the system and Bergen's technique for finding the detailed behavior of the 
limit cycles. Ina loop with a three-level relay, the width of the zero 
level may be adjiisted to eliminate the limit cycle. The response of 
the loop to inputs may be found by determining bounds on the con- 


Secutive numbers of pulses possible under certain conditions. 


The output sequences may be verified by a modified form of Bergen's 
method. For small inputs a modified step-by-step transient approach 
is best. This technique relies on detailed knowledge of the system's 
modes. 

3. Delta-modulated integrating gyros and accelerometers are 
useful inertial instruments. They deliver digital outputs in synchronism 
with a clock, and can be scaled as necessary to meet computational 
requirements. Delta-modulated instruments are fixed level torquing, 
which is much more accurate than analog torquing; and the measurement 
of integrated torque, which indicates the system output, is merely a 
pulse-counting process. 

4. The transformation of the A V vector from a body-fixed 
coordinate frame to an inertial coordinate frame is best accomplished 
through a direction cosine matrix transformation. The direct application 
of the matrix differential equation is most suitable for digital computation. 
Other methods considered for transforming vectors, while providing 
valuable insight into the problem, are computationally cumbersome. 

Oo. Since the instrument outputs are incremental in form, the 
digital differential analyzer is the most attractive choice for instrumenting 
the solution of the matrix differential equation. Their speed and compact- 
ness are attractive features for airborne computers. 

6. Integration rules with order higher than zero, using a larger 
increment, cannot be used to attain a slower computer speed with 
increased accuracy because of an irreducible uncertainty in gyro 
information. This uncertainty is a consequence of the digital form of 
the angular information. Thus, rectangular integration, being the 
simplest, is the best for this application. 

7. The transformation computer may be analyzed using frequency 
domain techniques. With the assumption that the angular rate about 
one axis is a rational multiple of that about another, the computer 
may be modeled as a sampled-data system, anda simple signal flow 


graph drawn for it. This signal flow graph provides great insight into 
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the computer operation, showing the dominant mode, and giving an esti- 


mate of truncation error. Quantization may be treated statistically. 
7.B Suggestions for Further Study and Development 


The areas of further study and development which can contribute 
to the realization of practical gimballess inertial reference equipment 


are listed below. 


1. Development of gyroscopes designed specifically for gimballess 
applications should be pursued. The severe angular velocity environment 
to which these instruments will be subjected provides enough design prob- 
lems to warrant a separate design. A torque generator providing high 
torque very accurately is the most important element in such a design. 

2. Moving-coil and flapper-type torque and signal generators should 
be developed. Their advantages of small size and greater sensitivity com- 
bine to make them an attractive choice for use in inertial instruments. 

3. The frequency domain analysis of the transformation computer 
should be extended to more complicated systems. A method should be 
developed which generalizes the results for rates which are not related to 
each other as rational multiples to cases where the rates are irrational 
multiples and to cases where the rates have varying relationships. 

4. The reversibility of the computation is an interesting point which 
deserves further study. The logic required in the computer, as well as 
the actual method of computation, should be considered. 

Do. stability in the large of the limit cycles of delta-modulated loops 


Should be investigated. 
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APPENDIX A 


DETAILED ANALYSIS OF A DELTA-MODULATED LOOP 
WITH A TWO-LEVEL RELAY 


The system to be analyzed is shown in Fig. A-l. This is the loop 
considered by many investigators with T. = 7 =1]second, D=1. The 
moding analysis has been carried out in Chapter II, Fig. 2-4. The 
Moding Plot indicates the existence of the 1-1, 2-2, and 3-3 modes. 
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Fig. A-l Delta-Modulated Loop to be Analyzed 
A.1 Zero Input Behavior 
A detailed analysis of the 3-3 mode follows: 
The output of the sampler is 
2 3 
+zZt - 
C(z) — ee (A -1) 


l1-z 


Here z is a unit delay, z = 5. Ge 


The output of the linear plant is 


A. TA, 
Aanz) = CQ) ae SGeTD) (A-2) 


The z-transform of A, giving the samples of A, is 


; A, vA, 
A(z) = -C(z) G(s) * l= a 











= = (1+z+z°) (1-29) (1-2) ace fey uf 
eo Set ee = 
1-z (1-2)? 1-z l-az (A-3) 
+ a + 1A. 1 i 
l-z l-z  £1l-az 
-T 
Where e s4c tire . 


If the 3-3 mode is possible, A(z) must be expressible in the form 


A +AztAgz +A w+ A yz +A 2° 
Aien= 2 Ml 3 2 
l- 2° (A -4) 





with the constraints 

A,790A,70A, 70 A,<0A,<0A, <0 (A -5) 
so that the postulated relay action is sustained. If Eq. (A-3) is to be put 
into the form of Eq. (A-4), the factor (l-@z) must be eliminated from the 
denominator of Eq. (A-3). For this factor to be eliminated, the numerator, 
M. (z), of the terms having (i-z°) (1-@z) as denominator must contain that 
factor. The numerator is 

M(z) = -(1-2°)*r41-2") TA, 


Hut 


a. (-2°)/ 12> + A. (1+z°)] (A-6) 
To 
Ss (1+A) eK Le ve z 
1+ A, 


If M(z) is to contain (1l-@z) as a factor ie must be such that 


or (A-7) 





This is one condition on the state of the linear plant for no transient in 


the limit cycle. With this condition satisfied, the output of the linear 
plant is 





A(z) Matiz”) eke op OTS ltzte +z °+2 *+2° (A is 
l-z (1-z) l-z l-z 2 2 
27 Wee ta-2*) (1-2) 
-A (A-8) 


l+a l-z 


The last term in Eq. (A-8) corresponds to M(z) after cancellation. 
Collecting powers of z in Eq. (A-8) gives 


Cae) ae A. 


3 
a - 207 + Soa T-T_+1)z 
q l+a@ l+a@ 


2 3 
+ (A - 2244 87-27 +1) 2° 











2 (ee Eos 
oa ince 3 
+ (A 5 ae 1) Z 
os l+a l+a 
207 ie 4 
+(A_ + 3 + 3 71-2T -7)z 
se l+a lta = 
Z 3 
+(A +e ty ==S,, 77)? 
Oita lta = (Aaah 


The numbers above ie == | seejmead ia 
(lez eA) = A. + (A, + .204)2 + (A,- 485 3)2- (A. - 1.195)z° 


+(A, - 1 304)2 7 + (A. - 3372" (A -10) 


Li 


By applying the inequalities Eq. (A-5) to Eq. (A-10) the following 
bounds on A, are found: 
.353 << A. < .837 (A -11) 


Thus, with initial conditions as specified in Eq. (A-7) and Eq. (A-11), the 
3-3 mode will occur without transient. 

Similar computatias lead to no transient initial conditions for the 
1-1 and 2-2 modes: 

l-l mode: 0< AG < .076 


A = a7 = . 463 (A-12) 


2-2 mode: 0<A< «415 


— = 762 (A -13) 





Further computations show that no other modes are possible. 
A.2 Transient Response 
A.2.a Impulse Response 

The impulse response of the loop is found by a modified step-by-step 
computation. After an input impulse, the output pulses remain one sign 
for a number of pulses determined by the strength of the impulse, as 
the increments of the output corresponding to the integrated impulse 
build up. As the output impulses pass through the zero order hold, 
they become a step. Thus, the appropriate calculation to find the state 
of the linear plant as the error dies out is a combination of step and 
impulse response. 

For anumerical example, assume that the loop is in the 1-1 mode 
centered at the origin. Let the input impulse be applied halfway between 
sampling instants, as the linear plant output is crossing from negative 
to positive. Let the imput impulse have strength S = 5.5. 

With t = 0 at the clock pulse just preceding the application of the 
impulse, the following equations describe the linear plant output 

A(t) = t - .038 -1.463(1-e') OS t 8.5 (A-14a) 
A(t)=t-.377+ 4.612 (te $5 St = (A -14b) 
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A(t) = -(t-1) + 2.403 + 3.80-e ty 1s ¢ (A -14c) 
Eq. (A-14c) changes sign between t = 7 andt = 8. Att = 8, the state of 
the system is 

A(8) = .792 and A(8) = -1 (A-15) 
These initial conditions lead to the 3-3 mode, as may be seen by 
comparison with Eq. (A-7) and Eq. (A-11). 

Note that the input integral is only determined to within three units. 
This is a quantization uncertainty due to moding. If the system had 
settled into the 1-1 mode, the quantization uncertainty would have been 
only one unit. 

A.2.b Step Response 

The first task in determining the step response of the loop is finding 
the upper and lower bounds on the number of consecutive pulses of one 
sign. To determine the upper bound on the number of pulses of one sign 
the maximum velocity and displacement of the linear plant output are 
assumed, and the number of pulses required to drive the plant across 
null is counted. If the magnitude of the input step is R = rD, the 
maximum rate that the system can attain is D(l+r). If a sampling instant 
occurs just before the system crosses null, the maximum possible 
displacement is D(l+r)T.. With these initial conditions on the linear 
plant and the linear plant being driven negatively by a step of magnitude 
D(-l+r), the time response of the system may be written as 

- 

A(nT_) = D(-l+r)nT, + Dl+r)T, + 2D7Q-e= 7 ) (A -16) 

If n is to be the upper bound on consecutive positive pulses, it must be 


the least integer such that 


A(nT _)<0 (A-17) 
For this to be the case, the relation between r, n, and the loop parameters 
must be 
me 
T -n— 
n-l-2 Te (l-e 7) (A -18) 


> 
—— OT 


Note that D does not appear in the above expression. For the numerical 


example above this relation is 


ee te alla 


10a: (A-19) 
Table A-1 is the result of substitution of various nin Eq. (A-19). 


tT > 


While this eliminates some possible sequences, the number 
remaining to be tested is still quite large; thus it is desirable to place 
a lower bound on the number of pulses of one sign which can follow a 
given number of pulses of the opposite sign. 

The computation to determine this lower bound is similar to that 
above. The initial conditions to be used are minima. The initial 
position is taken as zero plus, and the minimum velocity after n, 
negative pulses, resulting from maximum velocity at the previous 
Switching time, is 

T 


S 
= D(ltr)-2De 7 ™ 


A ,(n,) l (A-20) 


With the above initial conditions, the expression for the linear plant 


output is 


T 
n 


T 
= 
tT) (A -21) 


S 
—————— nN -_- 
A(nT.) = D(-ltr)nT, + 2Dr(I-e 7 ly (+e 


As before, if nis to be the required lower bound, it must be the least 


n for which Eq. (A-17) is valid. The required relationship is 





fl T 
ee eee x 
2D 7 (l-e 7 1)(l-e 7 °) (A -22) 
ae = 2 
n 


The ranges of r for various n and n, are shown for the numerical 


example above in Table A-2. 
The response of the loop for r = 5 will now be computed. Forr = = ; 


Table A-1l shows that while the feedback and input act together, the 


maximum number of pulses (which are negative) is two, and when they 


Ly 


TABLE A-l 


UPPER BOUND ON CONSECUTIVE POSITIVE PULSES 


range of r 








TABLE A-2 


LOWER BOUND ON CONSECUTIVE POSITIVE PULSES 
FOLLOWING A GIVEN NUMBER OF NEGATIVE PULSES 


aE EE AT 





range of r 


















































1 <r<, 2 “YT <re,095 | "1 <r<2.2 “| - 1 “<rera24 
2 i 22) ie a a0 =, 090sr—y2on a) —. 20 —neelo -.24 <r<, 15 
| 3 1) .456<r<.6 _ 20) SPO Noein Oe) lpn cress 
| 4 || .6 <r<.6ge .452<r<.579 | .4 <r<..533 | .38 <r<.518 
8 || .sgecr<.749 | .579<r<.856 | .533<r< .622_, 5 18<r<u6 11 


NOTE: Bounds on negative pulses may be found by changing the sign of r 


and entering the tables again. 


are in opposition the maximum number permitted is seven. For the 
feedback and input acting together, one is the lower bound on the 
number of negative pulses following any number of positive pulses, 
according to Table A-2. Howevey at least three positive pulses must 
follow one negative pulse, and at least four positive pulses must 
follow two negative pulses. 


Forr => the following is a partial list of sequences which have 


the required average: 


Period Sequences 
4T pe | 
Ss 
8T 41,254 Selle: 46;.2¢ 
12T. cs Me A is aie em) rg rors ra bets a rs le LL, 


Bide ds i Silat? ,.. Siglitieeugulis 2.9: 2s mene ele 
According to the bounds found above, the only sequences which are 
possible are 3,land6,2. The sequence . 3,1 will be verified ina 
manner Similar to that used in verifying modes above. 


The output of the loop is 


C(z) = arama. ee 2 (A -23) 


The output of the linear plant is 


* 


TA 


A 
O Q 
= - 4 ee 
A(s,z) = [R(s) - C(z)G,(s)] G{s) zi SiGe Te. 
The samples of the linear plant are given by 
~ le 7, = 
A(z) = D[ r-C(2) a-2)]] S$ ,.- 1 4-1 
(1-z) l-z l-az 
A+ 1A TA 
iano oF O 
l-z l-az (A -24b) 


If the postulated sequence is possible, then A(z2 must be expressible 


as 
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l-z 


A, +A zt Ayz” 4 Agz® 
—————————— (A -25) 


with 
(A-26) 


Am A, 7 0 A, ? 0A, <0 
For this to be the case, the factor (l-@z) must be eliminated from the 
denominator of Eq. (A-24). The numerator of terms with denominator 
4 
(l1-z°) (l-az) is 
[sea - Tee a ie TA w(l-z") 


a ae (tata -2° )] 


M(z) = 


T (1-z){(. 5-A.) 42424 








; 9 I. 5A, 2 
= - 7(.5+A.) (l-z) (l¢z+z~ - —————> z") 
-5tA) (A -27) 
For (l-@z) to be a factor of M(z), it is necessary that: 
ise 
te ongeta: 
.5+A 
O 
or 
ae Loe (l-apeeee) 
: Itata™ +a (A -28) 
With A , as in Eq. (A-28), A(z) may be written as 
4 2.3 Taz AG +TA, 
A(zy et —2Silez “We ete 22 02) a je *o 
1-z (1-2) 1-z l-z 
(A-29) 


- 41. (1- Zl lek a) Zac (lteta? ae ] 
Pere +a 


San COE powers of z yields 
(l-z*)A(z) = A. + (A +TA +. 57-. ST = — Fz 
S ltata“ta 


3 2 
HA + TA +57-T_- eee) a 
© S liteta“+a 
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AWl+at+e" ) 3 
)z 


lt+at+ta +E 


"id Ppetmome 


+ ee - - 
vy alee low = 1. oT. + (A -30) 
Substitution of the numbers of the example into Eq. (A-30) yields 
(1-2 ")A(z) = A +(A, +1. 312)z = (A+3.113)27 + A. - .25)2° — (A-31) 
Application of the inequalities Eq. (A-26) shows that for 3, 1 to be an 


output Sequence, the acceptable range of initial plant values is 
O>A,>.25 (A -32) 


Similar computations show that 6,2 is also a possible output 
Sequence, but that none of the others is possible. 

The above analysis, as noted in Chapter II, is useful only for 
inputs larger than about 10% of the saturating signal. For example, an 
input step of one-twentieth of the saturating signal would have an output 
sequence of some multiple of forty sampling periods, which implies a 
formidable computation. 

For small inputs, the modified step-by-step transient analysis 
described in Chapter II is the most useful. The key to the analysis is 
the fact that a small input may be viewed as a small disturbance of the 
prevailing limit cycle. 

For a numerical example, assume the system analyzed above 
is in the 1-l mode with zero average and that an input step of R= .01Dis 
applied at the sampling instant when the linear plant is positive. 

The plant is assumed to move under the influence of the input 
alone using the initial value of linear plant output and the velocity of 
the plant relative to the limit cycle as initial conditions. These initial 


conditions are 


A. = .038 Ao = 0 (A -33) 
and the equation describing the motion is 
_t 
A(t) = .O1Dt + .036 - .O01DT(l-e 7) (A -34) 


The mode will change when the plant output sample which is now 


negative becomes positive. When this occurs,t = (2K=1)T., and 


} 


Pall 


A[ @K+))T >.076, that is, the linear plant output sample gets outside 


its permitted range for sustaining the 1-1 mode [see Eq. (A-12)]}. 


Ss 
) 

(A -35) 
For K = @2, Al 2K+1)T J] = .07607, and the mode will change. The 


velocity of the(linear plant at this time is the drive velocity, .01, plus 


a 
A[ (2KH)T.] = . O1D(2K+1)T, + . 036 - olDt(1-e 7 (2h 


the limit cycle velocity, -. 463, [see Eq. (A-12)], and the position of 
the linear plant is .07607 - .076, since the plant output which was 
negative is now being considered. One sampling period later, the 


linear plant output is (shifting the time origin) 
A(T.) = -,.99 + .00007 + (-. 453 + . 99) ile) 
= See (A -36) 


The plant velocity at this time is 
A(T.) = -.792 (A-37) 
By comparing this state of the linear plant with the mode initial 
conditions, it is seen that the system will settle into the 3-3 mode. 
Again, the plant is assumed to move under the influence of the input 
alone, with initial conditions which are the actual plant position and 
the velocity of the plant relative to the limit cycle. Thus the initial 


conditions to be used are 


A, = -. Gol Ay = -, 792 - (-.905) =+ .113 (A -38) 
The equation of motion for the linear plant is (again shifting the time 
origin) 
a! 
A(t) = .01Dt - .651 + 7(.103)(1-e 7) (A-39) 


The mode will change when the linear plant output sample which is 


least negative becomes positive. This occurs att = (6Kt2)T.. At this 


time, the linear plant output will leave the region leading to the 3-3 
mode so that A[ (6K+2) T,] > - .353 [see Eq. (A-1])]. 


T 
-(6K+2)— 


A[ 6K+2)T_] = .01D(6K+2)T_ - .651 + 7(.103)(1-e io 
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For K = 3, A[6K+2)T.] = .348, and the state of the linear plant is 
A, = .005, and AG =+ .739, which leads to the 2-2 mode. This mode 
prevails for 80 periods, when the 3-3 mode appears. This mode 
persists for 50 periods, then 2-2 mode prevails for 50 periods. The 
average output has changed one unit in a hundred periods, giving the 
required average of .01. A plot of this is shown in Fig. A-3. 
A.2.C Effects of Imperfect Integration 

With imperfect integration in the linear plant, the integrator 
is replaced by a first order lag. In that case the linear plant in the 
forward path becomes 

G(s) ‘ (7)8 + ) (Tos + 1) (A -41) 

To show that a small input cannot be detected, a sequence with 
zero average and a small input will be assumed. The inequalities 
resulting from the output sequence will show the size of the input that 
cannot be detected. 


Assume an input R = rD and test the 1-1 mode. Then 





i l 
2 a et 
A(z) =Ir - i | = 1-75/7, _ 1-7, /T. 
l-z 1-z l-@)z l-a52 
A tT oAL AStTA, 
{ 1-7/7, r 1-7 JT 
l-a,Z 1-@.z (A -42) 
l 2 
T ap 
__s _ 5s 
where a =e Ty and a =e To 


If the postulated output sequence is to be sustained, A(z) must be of 
the form 


A +A,z 
Gl 


A(z) = 5 


1-z (A -43) 


with 
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A, >0 A, 0) (A -44) 


Again the (l-a@z) factors must be eliminated from the denominator of 
Eq. (A-42). The constraints on the initial conditions which lead to 


that elimination are 


l-a@ 





; I 
+ = 
A, Toh guests ita, 
; l-a, 
A, BE TAS = ee ee 
2 (A -45) 
With these conditions met, the output of the linear plant may be written 
as 
: l-a, i l-a@, 
9 llta, 2 l+a, 
(l-z-) A(z) =|r+ eee 
l 2 
til i | 
l lta, 2 +a, 
= 5) ee ae lz 
1 2 / (a-46) 


For the postulated mode to be sustained, the limits on r are 


: 1-a, , l-@, 
ll+a, ee) +a, 
l 2 


Similar calculations show that the magnitude of input which cannot 

be detected is that which would, in the steady state, drive the sample 
of the linear plant output closest to zero just to zero. For design, 

the dead zone for each mode must be found and the largest used as a 
limit. As was seen in the response of the loop for a small input, each 
mode was present, and the largest dead zone would have stopped the 


output. As the above discussion shows, elastic restraint causes an 
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error in the output. The average position of the linear plant and the 
elastic restraint constant combine to give an average force on the 
linear plant which is subtracted from the input. Thus, the output 
indicates the integral of the difference between the input and the 
average elastic force. 
A.3 Frequency Response 

The describing function of the loop shown in Fig. A-1l will 
be generated, following the method outlined in Chapter II. 

The input to the relay is assumed to be sinusoidal. In complex 


notation 

A(t) = Ae!” (A -48) 
Consequently 

c(t) = = el (A-49) 
The error signal is 


A. jot mo A j(wt-/G(jw) ) (A -50) 
E(t) = Gay ©  =—Tatjor ° psie), 


The feedback signal is 
oT 





sin S oe 
_ m4 2 (jwt- ) (A -51) 
B(t) = T wT = é 
Ss 
Z 


The input is the sum of the error and the feedback signal: 
R(T) = E(t) + B(t) 








oo wT 
A a sin —-=— Ss 
ere ae ee 
Ss 
2 (A -52) 


The describing function is the ratio of output to input amplitude 


=< 
1 


Cs te ne : 
om wT 


; , sin ee a 
a - JLGGe) + 4 ee Oo) 2 (A-53) 


lGGe)| T oa 


2 
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Fig. A-4 Low Frequency Describing Function of Loop 
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This function is frequency-independent at low frequencies if 
there is an integration in G(jw). A polar plot of this function is shown 
in fig. 2-4: 
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APPENDIX B 


ANALYSIS OF DELTA-MODULATED LOOP 


WITH THREE-LEVEL RELAY 


The same system as in Appendix A will be analyzed, with a three- 


level relay as shown in Fig. B-l. 





Fig. B-l Delta-Modulated Loop to be Analyzed 
B.1 Choice of Zero-Level Width 
As explained in Chapter II, the zero level of the relay should be 
chosen to be as narrow as possible while eliminating the limit cycle. 
There are two aspects of limit cycle elimination to be considered: the 


low-energy region, handled by a transient analysis, and the high-energy 


region, handled by a transient analysis, and the high-energy region, 
handled by a describing function analysis. Each analysis gives a lower 
limit on the zero level, and the larger is chosen as the size to be used. 
The low-energy analysis assumes that the linear plant output is 
at the edge of the zero level with zero derivative, and that one clock 
pulse occurs to drive the plant back to null. The zero level is chosen 
wide enough so that in an infinite time the other edge of the zero level 
is not reached. 
While the clock pulse is driving the plant toward null, the plant 


output is 
t 
INGE D, ted a lle (B-1) 
At t= T. 
. | Ts 
AT) = -D(l-e ~7) (B-2) 
After t = T, 
qT, t 
At) =" ~ a + At TD(l-e “T) - TD(l-e 7) (B-3) 
if A (co) > - A.: J must be chosen such that 
fev > ap Pap (B-4) 


This will be the result of any low-energy analysis, no matter what the 
dynamics, as long as the plant contains one integration. 
The high-energy analysis makes use of the describing function for 


the relay which is 





(B-5) 


This has a maximum, and hence = has aminimum atA = N24 | where 


NN2A) = TA (B-6) 


The size of the zero level is chosen so that there is no intersection 
of the - sy locus with the extended G(jw) locus. The magnitude of the 
G(jw) locus is critical for the 3-3 mode, the possible mode with the 
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largest amplitude. At that frequency 


G w w Da SG 
|. s Parc a 7 | (8-7) 
ig Otig 1)| LG op) 
© 1 
Equating |G(j —) land Tn yields 
A>tor 3 
e 2 


s \r ye 
a (3 7) (B-8) 
S | 


For slow systems (= <1), this width for the zero level is smaller 
than the low-energy requirement, and for faster systems it is somewhat 
larger. However, the factor multiplying DT, depends on the highest 
permissible mode. For the case being considered here ee = 7) Eq. 
(B-8) becomes 
A. > = DT. (2.08) = 1.04 (B-9) 

Since this is larger than the size specified by Eq. (B-4), it will be 
taken as the size of A. for this example. 
B.2 Transient Response 
B.2.a Impulse Response 

The impulse response is again handled in a step-by-step transient 
manner, recognizing that an impulse input changes the derivate 
of the input instantaneously by the strength of the input. 

For a numerical example, take the input impulse strength as 5.5 
and assume Zero initial conditions. At the end of one sampling period, 


the linear plant state is 


T 
8 
A(T.) = 5.5D7(l-e “7 ) = 3. 48 
T 
PT ils 5.5De" = = 202 (B-10) 


At this time the feedback signal appears; the linear plant output is 
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INDICATED INTEGRAL OF INPUT 
INTEGRAL OF INPUT 


Ts ITs lOTs 


Fig. B-2 Response of Loop to Impulse: R(t) = 5. 58(t) 
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DP 
a r 
A(nT.) = -DnT. + 5.5DT(l-e 7) + D7(5. sa * 1) 


T 
Ss 


(l-e ) (B-11) 
The feedback signal will cease when A(nT.) <A, orn = 6 in this 
example.. The state of the linear plant is then 


A(6T.) -9 


-l (B-12) 


A(6T.) 
The plant continues under the influence of the remaining velocity, but 


no feedback signal: 
oo 
A(t) = .5DT,=TD(1-e 7) 


and (B-13) 
A(*o) = -9DT_-TD =F 


Thus no negative pulse is generated. 

A plot of the loop's impulse response is shown in Fig. B-2. 
B2.B Step Response 

The step response of the loop with a three-level relay is computed 
in much the same way as that of the loop with a two-level relay. The 
bounds to be determined in this case are the number of consecutive 
output pulses of the sign of the input permitted and the number of 
consecutive pulses of sign opposite to the input. 

The first computation proceeds from the assumption that the 
linear plant output has travelled a full sampling period under the 
influence of the input alone at maximum velocity. This leads to initial 


conditions as the feedback signal appears 


A =A + DrT 
A =Dr 
O 
The linear plant output is then 
iv 
Aint) = -Dioanii + AW DrT_ Die T) (B-15) 


The maximum number of output pulses of the same sign as the input 
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is that n for which A(nT, A If this is the case, then 


T =T) qT. 
n- 7 (l-e 7?) 
> S 
2 nti (B-16) 


Substitution of various n into Eq. (B-16) yields Table B-1, which 
shows limits on r for a given maximum number of output pulses of 


the same sign as the input. 


Table B-l 
Limit:s on Size of Input,Given 
Number of Consecutive Output Pulses 


of Same Sign as Input 





The possibility of pulses of opposite sign from the input must be 
considered. If the feedback and input act in opposition for a full sampling 
period before switching occurs, the state of the linear plant is 

A = A_ - D(l-r)T 
O Cc Ss 
: (B-17) 
eM = D(l-r) 
Even for the smallest input, Ay > “AY so there will be no feedback 


signal initially. Thus the linear plant output is 


“<n — 5 
A(nT .) = DrnT oa A, - D(l-r)T. ~ WDil-e T ) (B-18) 
If ever A(nT .) < - Aw there will be a negative pulse; thus the expression 
if as ae 2A. 
ae . 7) =e 
iL. T. (l-e ) DT, 
= ntl (B-19) 
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is of interest. Clearly r has been assumed positive. For the present 
example, there can be no negative pulses for a positive input, and vice 
versa. 

The output for an input step R=5 D will now be computed. 
According to Table B-l, there may be one, two, or three positive pulses, 
but not four. However, the only possible output sequences are 1,1,0,0, 
and 2,2,0,0. These may be verified in the Same manner as the sequences 
tested in Appendix A. a 

For an input step of R = .01D, the only possible output sequence is 
1,99,0,0. That no more than one consecutive positive output pulse 
is possible is clear from Table B-1; thus, the above sequence alone is 
permitted. A short transient calculation shows that the first output 
pulse will occur at t = 1047; because of the width of the zero level. 
B.2.C Effects of Imperfect Integration 

As in the previous case, the effect of imperfect integration is to 
mask a smallinput. The limit on detectable inputs can be determined! 
by comparing the steady-state response of the linear plant with the zero- 
level width of the relay. Inputs which cause steady-state displacement 
of the linear plant less than A. will not be detected. Further, as before, 
an input corresponding to the average value of the linear plant output 
will be subtracted from the indicated output. 
B.3 Frequency Response 

As in Appendix A, a frequency dependent describing function may 
be calculated for the loop, beginning with an assumption about the relay 


input. Assume this is 





A(t) = Ael*t (B-20) 

The output of the loop is 

2 
city = 44/1 - (58 
T A (B-21) 
the error signal is 
= A ee ae re Ge) B-22) 
mere) ¢ 'GGen ° 


and the feedback signal is 
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B(t) = Bed) 





The input is the sum of the feedback and error signals, so that the 


describing function of the loop is 




















x 
4y/, | 4c 
T A 
= (A,) = wT =) Se Ss 
Soh A sin .~ e 2 
A_ _-j/G(je) 5 4,f- |< 2 
| G(Gjo)| 7 A oT. 
Z 
(B-24) 


This may be plotted as a function of R for each frequency. Unlike the 
describing function for the loop with a two-level relay, this is not freqency 
independent for low frequencies. 

Due to the zero level of the relay, there is a magnitude which R 
must attain, depending on frequency, before any output appears. It 


is required that 


A = RIG(u)I>A, (B-25) 
or 
Aa 
R° TGGo)l (B-26) 


for an output to occur. 
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APPENDIX C 


ANALYSIS OF TRANSFORMATION COMPUTER 


C.1 Roots of Characteristic Equation 


The roots of the computer's characteristic equation for two dif- 


ferent inputs will be determined using the root-locus method. 
C.l.a Equal Rates on All Axes 


In this case W) = Wo =W9, SO Z) = Zo = Ze = Z. Eq. (5-27) becomes 


a Er. ee (C-1) 


(= 26 (1 - 2)? 


Weer yan? 


As is shown in Fig. C-1l the roots of this equation remain on the 
unit circle. This indicates a non-decaying oscillatory behavior, as 


expected of direction cosines. 
C.l.b Different Rates on Different Axes 
Suppose that the angular velocities are W) = 2Wo, We = QO. The 
difference equations are then 
C,,(2n) = C,,(2n-1) ~ AGC, 2(2n-1) 
C,,(2n+1) = C,,(2n) 
C,5{n) = C.9(n-1) + AGC, .(n-1) 
C,,(2n) = C,3(2n-1) a A@C,,(2n) = AGC. 4(2n) 


C,,(2n4) = C,,(2n) - AC; (2n+1) (C-2) 





AGW AVY AY S74 = 


(1-Z)2 ’ 


Age'= ./3A@ 


Fig. C-l Root Locus For Equal Rates on all Axes 
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To transform these equations into the frequency domain, it is 
OO 


necessary to compute y C.4(2n)2”, which is the transformed sequence 


n=0 
of every other value of C5. This is the sequence which is fed into the 
C,, node (register). 


2nt+l 


on n 
2C,3(2n)z 2C,,(n)z 2C.(2nt1)z 


2n+l 2n+l 


n 
ZC, (n)z - 2~C,(4n)z a AGC, .(2n+1)z 


2n+l 


n . an 2 2n+l 
2C,,(n)z - z2C.,(2n)z + A82C,.(2n)z + AO 2C.,(2n)z 


EC, a(n)z™ - 2(1 - A@*)EC, 4(2n)z = A@zZC..(2n)z" (C3) 


ae 
ZC, 9(2n)z°" = EC, o(n)z™ - ZC, 9(2ntl)z"* 























2n+l 2n+l 
= EC,9(n)z" - ZC, (2n)z nT" = A®'EC, 4(2n)z" 
74 | | n . 2n 
Be (an) z = =C.,(n)z - A@z2C..,(2n)z 5 (C-4) 
i2 i2 13 
l+z 
EC, 4(2n)2™ = ZC, 4(n)z” - a (1 + Ag ee I] Bc, (onic! pos Cuan 
to nye 1) i2 
l+z l+z 
a (C-5) 
i + 2[1+A0" (1 + =e} 
l+z 
In z-transform theory, the z-transform of a function is defined as 
CO OY 
he n _ on 
C,.(2) =) en (Ref. 61, Ch. 1). Let C32) » Cagis Me . Then 
n=0 n=0 


Eq. (C-5) may be written 
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Az 








| C53(2) = : C59 (2) 
C:3 Z) = (C-6) 
fi +2}[1+a0° [t+] 


From the difference equations are obtained the.relations 











Cc! : 
me — C,,(z) (C-7) 
C.(z) = - =< C!a(z) (C-8) 
A8z 
Cag) se C, (2) (C-9) 


AQ 


=a s 





Gra) “—S>- Cpl) - Ch (z) (C-10) 


Combining Eqs. (C-6) through (C-10) gives the characteristic 


equation for the computer: 








222 
ae? (1-202 
2 2 
1 + A9 z ee oo) ee er; (C-11) 
2 2 2 Z, 
(1 - z) (l-z)Q+z)j/l+¢z{1L+A0° [1 + 
[+2 fh +a0"(1 +f] 


Fig. C-2 shows the root-locus diagrams used to find the 
singularities of Eq. (C-ll). There are four important zeros of this 
equation, marked with squares in Fig. C-2(b). The pair of roots near 
unity represent the dominant mode of the system. The pair near -l 
are oscillatory with magnitude 1 ae . This gives rise to a time 
sequence (1 - ey" Since the magnitude of the mode is O(A®@) the pair 
of roots gives rise to an error growing as nae’, This leads to an 


increment of error in about one revolution. 
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SINGULARITIES FROM 


T 1 AND —| ; 
al ARE 0(A@2), 
ReZ 
{- 2 
| — ya 





ee a SS 
(14Z (142 [14 002(1+ 4 ) J) 





Ag2z% 
* es ( ‘ Sez : 
(1-z)2 (1+Z)(1+Z[1+A0*[I+4, ]) 


Fig. C-2 Exact Root Locus for W1 F205, 0,=0 
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C.l.c Approximate Analysis of Multirate System 


For W) = 2Wo, Wa = 0, Eq. (5-22) shows k = 1, @ = 2, and iq. (5-23) 
shows Z) = Zo = hae Substitution of these relations into Eq. (5-25) 
yields Eq. (5-21). The characteristic equation corresponding to this 


setis 


2 2 
A’ z AS z 





2 
(ezy? Bac 


Comparison of this equation with Eq. (C-1L) shows that it is approxi- 
mately correct. The root locus diagram for this equation, shown in 
Fig. C-3, is practically the same as for Eq. (C-ll). The extra 
dipoles appearing in Fig. C-2 have negligible effect. 


C.2 Trumeation Error 


Following the method outlined in Chapter V, Sec. 5.D.3, an 
analysis of the truncation error made by the computer will. be made, 
first for equal rates on all axes, then for W) = 2Wo; We = 0. In each 


case, the initial condition C(0) = I, will be used. 
C.2.a Equal Rates on All Axes 


The exact solution to Eq. (5-1) in this case is 





2 
lis 44 . Sones Sal 
3 Z 3W eo oy 3 
Ww 
2 
C(s) = aA 18 ane 2 ee mene (C-13) 
Ss SHO 31.2 3w «Ss 33 
a W : 
3W 
= {Be a, 
3w «Ss 8 3H 8 3 oe 


The time solution is then 
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Fig. C-3 Approximate Root Locus for oO 2u =0 
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1+ 2 cos V3 wt 1 - cos V3 wt - V3 sin V3 ut 1 - cos V3 ut + V3 sin V3 wt 
cit) «4 1 - coe V3 ut + V3 gin wt 1+ 2 cos V3 wt 1 - cos N3 wt - V3 sin V3 ut (C-14) 


1 - cos V3 wt - V3 sin wt 1 - cos V3 ut + V3 sin V3 ut 1+ 2 cos V3 wt 


According to the flow graph shown in Fig. 5-4, the computer 


solution, without quantization is 



































A672 a6 A6z a6 A062 
1+ -—— + + 
(1-2)? ~ J-z (1-2)? l-z (1-2)? 
Toa) Adz A672 b672 Aé 2 
C(z) <= + 5 1+ 5 = a (C-15) 
1 + 340% + 46%)z2| 1-2 (1-2) (1-z) l-z ~ (I-z) 
(1-z) 
A6z 6°22 d6z d6%2 A672 
: + 5 7 1+ 5 
l-z = (l-z) l-z (l-z) (l-z) 


The normal modes of this determinant have unit magnitude and 


argument 


ry 


nf + Se - Sell ig 22° 


AQ’ Bgne 2 ae 2 2 et (C-16) 


$3405 "ge 


| 222 eee 


2 4 


Thus the sequence corresponding to Eq. (C-15) is approximately 





1+ 2 cos n S80! 1 - cos n 4@' - VJ gin n A@' 1 - coe n 46' +H sin n 48! 
a@ 
1 - 22 = 
C(n) = 3 1 - cos n 48 + V3 sin n J8' 1+ 2 coa n 486! 1 - coe n A6@' - sin n 48! Ce 1) 
3 
1 - cos n 4@' - V3 ain n 48' 1 - cos n 46' + 3 einn 46! 1+2 cos n 46' 
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The truncation error arises in this case from the fact that the 
argument of the roots near +l is not quite correct. This leads to an 


error which grows as nde’, 
C.2.b Multirate Systems 


With W) = 2Wo» Wa = 0 and C(0) = I, the exact solution to 
Eq. (5-1) is 





s? + we n 48 

5 2 ; 

4 

2 
so +8 
C(s) - l/s __ .4 a a= (C-18) 
E 5 2 i 
+ 
5 2 
—W 
-.45 3S l 
W W 
The corresponding time solution is 
Ls (4 + og wt) ss (1 - cos BE wt) oye sin MEE 
5 2 5 2 N5 2 
C(t) 2 - Bee st) ti1+4cos N35 Oty = pti NB ut 
5 2 5 2 N5 2 
1. N5 20) 5: V5 
— sin — wt - — sin — wt cos —— wt 
N5 2 N5 y) 2 
(C-19) 


Under the same conditions the model of the computer is as 


shown in Fig. C-4. 


147 





Fig. C-4 Computer Model for WF 2W5, 0270; C(0)=I 
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2 2 
(1-2)? (1-z)(1-z7) 1-22 
Cr A6*2 A67z A6 
oe nae re: eae 1 ale (C-20) 
i+ AG"z, 46"z (l-z)(1-z") (1-2) (l-z) 
(1-2)? (1-27) 
_ Adz Az 





l-z l-z 


There are essentially two parts to the time solution of 
Eq. (C-20). The first, desired part of the solution arises from 
the roots near +l. The roots near -l give rise to extraneous, un- 
stable terms. The roots near +l have magnitude 1 te EE) and 
argument A@'= BLEW) The time sequence represented by Eq. 
(C-20) is approximately 


A (4 + cos HAG!) f (1 - cos n4é@') = sin nAQg' 
5 5 N5 
C(n) = len SY ay - cos nA8') (i + 4 sos nA@') 2 sinnA@! 
g 5) 5 N5 
1 ssinnAe') 2 sinnAe’ cos nAé@! 
N5 N5 
Z 
= 0 i 0 Ag” 0 \ 
4 2 8 \ 
2 l ae” — ae? 
+nA@ 0 0 — |-| —— — 0 (C-21) 
2 8 8 
L L 1B 0 0 0 
2 2 4 


Again, an estimate of the truncation error may be had by com- 


paring Eq. (C-21) and Eq. (C-14). The error arises from the presence 
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of the poles near -l, and the fact that the roots near +l are not on the 
unit circle. The truncation error grows as nde? in each direction 
cosine, It should be noted that the approximations made in the analysis 


precluda:taking Eq. (C-21) as the exact computer solution. 
C.3 Quantization Error 


The power spectral density of quantization may be calculated 
following the procedure indicated in Chapter V, Sec. 5.D.4. However, 
perhaps a better indication of the error due to quantization is the mean- 
squared signal appearing at any node due to quantization noise. This 
will be done for the case of Wis = 2Wo; We = 0, using both hysteretic and 


non-hysteretic quantizers. 
C.3.a Non-Hysteretic Quantization 


The model to be analyzed is shown in Fig. C-5. The expression 


for power spectral density at each node.is: 





y, 2 2 
1 +49 Z 4 AQ z : AQz 
®.. = (1 - ae (1 - ae -z) 1l- z° Ae* 
a, 
y) y, 12 
1 + 48 Z "7 AQ z 
(l-z)* (1-27)? 
2 2 
he Adz fi As’ z 4 Oz 
242 2 l1-z 2 
e 2) (a2 (1 - 2b a, AS (C-23) 
* Ae- ig 7 12 
| -+ -S25 2. eee 
(1 - z)° (1 - Zo 
ie A®Oz 7 A8 2 
1g2 1-2 ee 
33” 2 2 12 
A As z a AQ-z 
(eres C= Ze 
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Fig. C-5 Computer Model with Non-Hysteretic Quantization 





Fig. C-6 Computer Model with Hysteretic Quantization 


hey. 


Since the magnitudes ofthe roots are very close to unity, the 
mean square value will be very large. Thus, only n terms in the sum 
of squares will be considered. Note that the contributions from the 


roots near +l are oscillatory, so that the mean square may be taken 
2 
th 1) 
quantization noise at the ij node. Carrying through the computations 


as half the amplitude squared. Let qs. be the mean squared value of 


indicated above gives for mean squared value after n steps 


ae SS 


qé(n) = 096 nae? 


ao(n) =1,2 Oi 


(C-24) 


a’ (n) = 0188 nAe* 


The rms quantization noise at any node, then, is on the order of 


mae 


C.3.b Hysteretic Quantization 


The model of the computer with hysteretic quantization is shown 
in Fig. C-6. Note the power spectral density of the noise is four 
times as great as the non-hysteretic noise, due to the larger quanti- 
zation box, 

The self loops may be absorbed into the other paths, giving a 
tae tome tl “> z) multiplying each path gain. The expressions cor- 


responding to Eq. (C-23) are: 


2 
2 2 
ae Az ” AG z _ Lez 
(= 2)*(L- 42)° ot SY = edly ez) -z*)d-=2} 
; 2 2 2 | ae 
® = = 
ul 
Neen Ae~z 2 
eee ae 


(Wee) = eee 2 yee 
Zz 2 
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ee Adz i AQ'z a 
Q=27)70- ley @-z7%)a-2)0-L2) a-zja-!z) 


Az 









7 2 Q 9 Ae" 
ae | errr paid fe 
¥ Ae“z Az 3 

a a 
eee 777i neaai= 2 2)" 
2° 2 
2 
pic A8z _ Ad 
Seo 2) a-za-z) 
_ - 2 
ae g g a (C-25) 
m A.8°z A0-z 3 
I sp - ot 
= ay" Oa eee Olea 


2 a 


Performing a sum-of-squares. computation similar to the above 
gives: 


a(n) = .384 nA oh pela 8) Ae* 


q’y(n) = 4.8 nA@* + 5.23 Ae” 


q’(n) = .0752 nA@~ + .158 Ae” (C-26) 


The rms error is a factor of two higher than before, anda 
constant has been added. 
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A. Delta-Modulated Loops 
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aofit 1946 a Paris, République Francaise, Ministre de l'Industrie et 
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du Commerce, Service de la Propriété Industrielle. 


French patent of delta modulation with circuit 
diagrams and waveforms. 


7. Robert Hayum, Compensation of a Digital Integrating Accelerometer, 
M.S. Thesis, M.I.T. (MIT/IL Rpt. T-209A) June, 1959. 


Compensation is accomplished by using torque 
generator microsyn as torquer and sensor. 


8. Patrick B. Hutchings, Compensation of .a Sampled Data Contactor 
Servo, M.S. Thesis, M.I.T. (MIT/IL Rpt. T-222) May, 1959. 


First attempt to handle inputs. Heurtstic 
discussion provides good insights. 


9. K. Izawa and L.A. Weaver, 'Relay-type Feedback Control Systems 
with Dead Time and Sampling,'' Trans AIEEE, Part II, LXXVIII (1959), 


49-53. 


Phase plane analysis of delta-modulated loops 
with dead time. Attempts to determine modes. 


10. Rudolph E. Kalman, Phase Plane Analysis of Non-Linear Sampled 


Data Systems, M.S. Thesis, M.I.T., 1954. 


Essentially a computer study of a delta- 
modulated loop, with emphasis on stability 
in the large. 


ll. F.J. Mullin and E.I. Jury, "A Phase Plane Approach to Relay 
Sampled Data Feedback Systems,'' Trans AIEFE Part II, LXXVI111 (1959), 


517+ 52 3 es 
Phase plane analysis with a nice method for 
reintroducing time into computation. 


12. J.J. Rocchio, Jr., Analysis of a Digital Integrating Pendulous 
Accelerometer, M.S. Thesis, M.I.T. (MIT/IL Rpt. T-176) June, 1958. 
13. Frederick A. Russell, Design Criterion for Stability of Sampled 
Data On-Off Servomechanisms, Ph.D. Thesis, Columbia University, 


June, 1953. 


Like Chow (Ref. 3), Russell introduces variable 
delay due to sampling. Uses Bode plots. 


14. J.F. Schouten, F. de Jager, and J.A. Greifkes, 'Delta Modulation, 
A New Communications System for Telecommunications' Phillips Tech. 
Rev., ile 952), 2 meco0- 
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Discussion of delta modulation and its 
application to transmission. 


[S.seocuert C., Schwartz, Analysis of sampled Data Relay Servomechanisms 
with Zero and Step Inputs, M.S. Thesis, M.I.T. (MIT/IL Rpt. T-284), 
May, 1961. 


First use of Bergen's method to verify modes 
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technique to verify step input behavior. 


16. Joseph J. Sliwkowski, A Study of a Delta Modulator, M.S. Thesis, 
M.I.T. (MIT/IL Rpt. T-287), May, 1961. 
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Describing function for amplitude only derived. 
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limit cycles not eliminated. Discussion 
of logical compensation. 


B. Inertial Instruments 
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The Floating Integrating Gyro and its Application to Geometrical 
Stabilization Problems on Moving Bases, A. Sherman F’. Fairchild 
Fund Paper (FF-13)(Institute of the Aeronautical Sciences, New York), 
he See 


General discussion of gyros and their application 
to platform stabilization. Gyro performance is 
considered in some detail. 
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B. Inertial Instruments 

21. Charles S. Draper, Walter Wrigley, and Lester R. Grohe, 

The Floating Integrating Gyro and its Application to Geometrical 
Stabilization Problems on Moving Bases, A. Sherman F, Fairchild 
Fund Paper (FF -13)(Institute of the Aeronautical Sciences, New York), 
he) ay 


General discussion of gyros and their application 
to platform stabilization. Gyro performance is 
considered in some detail. 
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22. William E. Fellows, "Vibration Effects on Gyroscopes and 
Accelerometers, '' Minneapolis Honeywell Regulator Co. Aero 
Document R-Ed 29004-1A, November, 1958. 


Detailed discussion of many errors caused 
by vibrations. 


23. Richard H. Frazier, Analysis of Four Pole Torquer with Full 
Normal to Gap Faces, MIT/IL Rpt. E-700, March, 1958. 

24. Richard H. Frazier, Four Pole Alterating Current Signal Generator 
with Flapper Type Rotor, MIT/IL Rpt. E-703, April, 1958. 


Each reportis a brief discussion, with numerical 
examples, of theory of flapper-rotor instruments 
in different applications. 





20. Albert P. Freeman and Peter J. Palmer, Abridged Nominal 
Values and Specification for 2FBG-2 Series Floated Gyro Unit, MIT/IL 
Doc. Gi-INU-Uzzr BG-=2C, duly 1959: 

26. Philip J. Gilinson Jr., William Denhard, Richard H. Frazier, 

A Magnetic Support for Floated Inertial Instruments, A Sherman, 

F,. Fairchild Fund Paper (FF -27)(Institute of the Aeronautical Sciences# 
New York), May, 1960. ’ 


Discussion of theory of magnetic suspension 
system used in many inertial instruments. 
Numerical examples provided. 


27. William B. Haff and Melvin Meltzer, Effects of Angular Vibration 
on the Performance of a Single-Degree-of-Freedom Integrating Gyroscope, ; 
M.S. Thesis, M.I. T. (MIT/IL Rpt. T-257) May, 1960. 


Theoretical analysis of vibration effects on 
gyros. Coning and sculling discussed. Series 
solution for equation of gyro motion proposed. 


28. Robert K. Mueller, Microsyn Electromagnetic Components 
(Instrumentation Laboratory, M.I. T.),December, 1952. 


Fundamental discussion of microsyns in 
their various forms by the inventor. 


29. Herbert Weinstock, A Study of the Response of the Single-Degree- 
of -Freedom Integrating Gyroscope to Angular Vibrations, MIT/IL 
Rpt. E-885, January, 1960. 


Consideration of simultaneous angular 
vibrations about several instrument 
axes at once. 
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C. Coordinate Transformations 
30. Herbert Goldstein, Classical Mechanics (Cambridge, Mass., 1951). 


Chapter IV has a complete discussion of matrix 
transformation, Euler angles, and Cayley-Klein 
parameters. Good annotated bibliography. 


31. Sir William R. Hamilton, Elements of Quaternions (Dublin, 1856). 


This treatise introduces quaternions as an 
extension of complex numbers. 


32. A.S. Hardy, Elements of Quaternions (Boston, 1881). 


Basic text on quaternions with geometrical 
applications. 


33. H.W.L. Hime, The Outlines of Quaternions (London, 1894). 


Quaternions introduced as the quotient of 
two vectors. Various forms given. 


34. Sir Harold Jeffreys and Bertha Swirles (Lady Jeffreys), Methods 


of Mathematical Physics. (Cambridge, England, 1956). 


Chapters 3 and 4 cover matrices and tensors very 
well. Each chapter is headed by an apt quotation. 


35. Myron Kayton, Coordinate Frames in Inertial Navigation, Ph.D. 
Thesis, M:I. T. (MIT/IL Rpt. T-260) August, 1960. 


General discussion of coordinate transformation 
theory in Chapter 4 using tensors. 


36. P. Kelland and P.G. Tait, Introduction to Quaternions (London, 1873). 


Written ''for those who desire to become 
mathematicians. "' 


37. J. Halcombe Laning, Jr., The Vector Analysis of Finite Rotations 
and Angles, MIT/IL Rpt. 6398-S-3, September, 1949. 


An interesting development of algebra required 
for the use of vectors to represent angles and 
rotations. 


38. Philip M. Morse and Herman Feshbach, Methods of Theoretical 
Physics, 2 vols. (New York, 1953). 


Chapter 1 develops some transformation properties 
and contains a brief mention of quaternions. 


39. William R. Weems, An Introduction to the Study of Gyroscopic 
Instruments. (Department of Aeronautical Engineering, M.I.T. )5 


January, 1946. 
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Section 8 contains a nice discussion of Euler angles. 
40. E.T. Whittaker, A Treatise on the Analytical Dynamics of 
Particles and Rigid Bodies, qth ed. (Cambridge, England, 1959). 
Se — 


Chapter 1 contains discussion of rotational 
transformations using quaternions, Cayley- 
Klein parameters, and Euler angles. 


D. Numerical Analysis and Computers 
41. Anonymous, ‘Modern Computation Methods; } Notes on Applied 
Science, National Physical Laboratory (of England)(New York, 1957). 


Complete discussion of numerical analysis as 
applied to computers. Excellent extensive 
bibliography. 


42. Vannevar Bush, ''The Differential Analyzer," J. Franklin Inst., 
CCXII (1931), 447-488. 


Detailed physical description of one of the 
first analog computers. 


43. Vincent Rocco DeMarco, The Effect of Increment Representation 
on the Accuracy of the Digital Differential Analyzer, M.S. Thesis, 
M.I.T. (MIT/IL Rpt. T-158), May, 1958. 


Complete discussion of errors in DDA with 
various increment representations. Analysis 
in time domain. | 


44, George F. Forbes, The Digital Differential Analyzer, gee 
ed. (Pacoima, Calif., 1956). 

Discussion and analysis of DDA. 
45. Laurence J. Gitten, Incremental Computer Solution of a Coordinate 
Conversion Problem, M.S. Thesis, M.I.T., 1959. 


Study of a coordinate transformation using 
a general incremental computer. 


d 
46. David R. Hartree, Numerical Analysis, 2’ -ed.(Oxford, 1961). 


A standard reference, includes a section on 
organization for automatic computation. 


47. Francis B. Hildebrand, Introduction to Numerical Analysis 
(New York, 1956). 


Standard text on numerical analysis. 
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E. Miscellaneous 
96. Lewis Carrol, Through the Looking Glass, and What Alice Found 
There, illus. John Tenniel (New York, 1941). 


The quote is in Chapter I. The explanation, 
by an early astronaut who had trouble on the 
launching pad, is in ChapterVI. 


57. Harold Chestnut and Robert W. Mayer, Servomechanisms and 
; d 
Regulating System Design, 2 vols. Vol I, 2°" ed. (1959), Vol. II (1955). 


Comprehensive text on control system theory 
and practice. 


58. Ralph J. Kochenburger,''A Frequency Response Method for 
Analyzing and Synthesizing Contactor Servomechanisms,'' Trans AIEE 
Part I, LX1 (1950), 270-284. 


The describing function presented and applied 
by its originator. 


59. Albrecht L. Kosmala, Feasibility Study of a Gimballess Inertial 
Space Reference, MIT/IL Rpt. R-274, April, 1960. 


Early attempts at analysis of gimballess 
inertial reference equipment. 


60. Samuel J. Mason and Henry J. Zimmerman, Electronic Circuits, 
Signals, and Systems (New York, 1960). 


Chapter IV contains detailed discussion 
of flow graph analysis. 


61. Robert W. Sittler et al., Pulsed Data Systems (Course Notes for 
6.54, M.I.T., Spring, 1960). 


Course notes for sampled data systems. 
Includes material on sampling, signal 
recovery, compensation and statistical 
theory of quantization. z =e °-s 


62. Julius J. Tou, Digital and Sampled Data Control Systems 
(New York, 1959). 


Sampled data text. Good section on compensation, 
especially with digital computers. z =e Ss. 


63. John G. Truxal,Automatic Feedback Control System Synthesis 


(New York 1955). 


Standard text in the field covering practically 
all aspects of automatic control. 
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48. Alton S. Householder, Principles of Numerical Analysis, 
(New York, 1953). 


Numerical methods for sets of linear equations, 
differentiation and integration. Excellent 
bibliography and bibliographical notes. 


49. Harold E. Maurer, Error Analysis of a Digital Differential 
Analyzer, E.E. Thesis, M.I.T. (MIT/IL Rpt. T-129),May, 1957. 


Attempt to use Widrow's (Ref. 55) statistical 
characterization of quantization. 


90. Harold E. Maurer, An Approximate Analysis of Error Propagation 

in a Digital Differential Analyzer (MIT/IL Rpt. E-698), March, 1958. 
Extension of Ref. 49. 

51. John M. Salzer, Treatment of Digital Control Systems and 

Numerical Processes in the Frequency Domain, Sc.D. Thesis, 

M.I.T,, August, 195l. 


A clean way to handle numerical analysis. 

Presents a design method for computers. 

52. James B. Scarborough, Numerical Mathematical Analysis, ue 
ed. (Baltimore, 1958). 


Standard work on numerical analysis. Full 
of numerical examples. 


53. Wallace E. Vander Velde, ‘Computers in Inertial Navigation 
Systems'"', Section V in Inertial Guidance - Terrestrial and 
Interplanetary, (Course Notes for 16.45S, M.I.T., Cambridge, Mass., 
Summer, 1960). 


A good basic discussion of all three types 
of computers, their operative advantages 
and disadvantages. 


54. E.T. Whittaker and G. Robinson, The Calculus of Observations, 
4‘) 44. (Glasgow, 1958). 


Covers most numerical analysis topics as well 
as curve fitting and Fourier series. 


595. Bernard Widrow, A Study of Rough Amplitude Quantization by 
Means of Nyquist Sampling Theory, Sc.D. Thesis, M.I.T., June, 1958. 


A detailed analysis of quantization through 
a Statistical model. 


/ 
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64. John G. Truxal, ed., Control Engineer's Handbook (New York, 
1958). 


As the title implies, this is written for 
the practicing engineer. Comprehensive 
coverage. 


65. Albert D. Wheelon, ''Free Flight of A Ballistic Missile, 
ARS J., XXIX (1959), 915-926. 


One of the few places that range derivatives 
are discussed in the open literature. 


66. Henry F. Blazek, ''The Performance of Inertial Components 
on an Unstabilized Base", Presented at ARS Controllable Satelites 
Conference, M.I.T., 30 April 1959. 


A brief listing of the effects of unstabilized 
bases on instruments. 
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